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High  Power  Josephson  Effect  Sources 


Final  report  on  contract  ^  F1962890K0033 

The  goals  of  this  contract,  were  the  development  and  demonstration  of  the 
terhnokigies  needed  for  snbmillimeter  wave  sources,  based  on  arrays  of  Josephson 
jun'  fKuis.  with  power  levels  above  one  niilliwatt.  .As  a  demonstration  of  this  technology 
It  was  planned  to  build  a  source  with  about  100 /i\V  of  power.  These  goals  required 
several  <idvances  in  labrication  te<  hnology  and  design  concepts. 

Hiilii  power  retpiires  that  rhi‘  /,  of  the  junctions  have  a  value  of  tens  of 
iiitihamps.  .Achieving  this  wiiih'  keeping  the  juintion  size  small  enough  to  avoid  flux 
flow  instabilities  re<pures  values  of  ./,.  approaching  lOOk.A/cm*.  Near  the  start  of  this 
contra<'t  we  demonstrated  rh«'  ability  to  fabricate  .\b/A10j./Nb  tnlavers  having  this 
range*  of  ./,.s  ami  to  use  these  trilayers  as  part  of  a  cotnplete  cncuit  fabrication 
fechuology  to  make  sources  using  parallel  lua.sed  junctions  having  jjower  levels  in  the 
l/j\V  range.  This  work  w;is  presented  at  the*  1990  .Apiilied  Superconductivity  conference 
Ij  and  IS  described  in  detail  in  appendix  .A. 

The  array  dc.'sigii  usc'd  in  our  sources  iucc-rporates  the  junctions  in  a  long  niobium 
niicrostrip.  .At  frequencies  approaching  the  gap.  the  losses  in  this  inicrostrip  become 
significant.  .An  accurate  method  of  measuring  these  losses  in  the  submillirneter  wave 
range  is  needed  to  measure  the  loss  in  N’b  and  other  potential  stripline  material  as  part 
of  materials  fabrication  development.  .Al.'-o  an  accurate  measure  of  the  loss  is  needed  in 
assessing  the  operation  of  the  array  source  in  order  to  compare  the  measured  power 
from  the  source  with  that  predicted.  Such  a  technique  WcLs  developed  as  part  of  this 
contract.  Here  single  junctions  were  incorpo’-ated  in  series  with  a  niicrostrip  resonator 
made  of  the  material  whose  loss  v/as  to  be  measured.  The  junction,  which  was 
resistively  shunted,  acted  as  both  the  generator  and  detector  of  the  submillimeter  wave 
radiation  in  the  microstrip.  The  amplitude  and  frequency  of  the  resonance  can  be 
determined  from  the  perturbation  to  the  junction's  I-'V’  curve.  This  gives  both  the  real 
and  imaginary  parts  of  the  inicrostrip  impedance.  These  results  were  published  in  1991 
[2]  and  appear  in  detail  in  appendix  B. 

In  work  under  previous  contracts  and  that  of  Ref.  [1],  junctions  of  the  oscillator 
array  were  placed  as  wavelength  intervals  along  the  microstrip.  This  ‘quasi-lumped” 
de.sign  preserved  the  correct  phase  relationship  among  the  junctions  to  ensure  stable  in- 
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phase  locking.  However,  the  large  amount  of  microstrip  required  would  be  both  too 
lossy  and  take  up  too  much  space  for  higher  power  arrays.  As  a  result  of  simulations,  it 
was  found  that  groups  of  junctions  covering  about  1/4  wavelength  could  be  placed  at 
each  half  wavelength  interval.  A  small  “distributed”  axray  based  on  this  design  was 
fabricated  and  tested.  It  showed  fully  coherent  operation,  demonstrating  that  the 
concept  worked  as  predicted.  These  results  were  presented  at  the  1992  Applied 
Superconductivity  Conference  [3]  and  are  described  in  detail  in  appendix  C. 

An  important  feature  of  any  source  is  its  linewidth.  Theoretical  analysis 
indicated  that  the  contemplated  high  power  subrnillimeter  wave  sources  would  have 
linewidths  of  around  10  kHz.  which  is  sufficiently  narrow  for  a  wide  range  of 
applications.  These  theories  had.  however,  only  been  tested  in  lumped  arrays  in  the 
microwave  range  since  the  submillimerer  wave  .sources  used  on-chip  detectors.  Work 
under  this  contract  was  therefore  carried  out  to  directly  measure  the  linewidth  of  the 
array  described  above  (App.  C).  Two  distributed  arrays  of  10  junctions  each,  along 
with  an  SIS  mixer,  were  fabricated  on  a  single  chip  forming  a  combined  source-receiver. 
The  mixing  product,  with  an  IF  of  about  10 GHz.  was  coupled  out  of  the  cryostat  to  a 
spectrum  analyzer.  The  linewidth  was  measured  as  a  function  of  temperature  from  1.5 
to  4.2  K  and  fc>und  to  be  within  5091  of  the  theoretical  prediction.  This  result,  which 
was  published  in  1993  [4l  an<l  is  described  in  detail  in  appendix  D.  gives  strong  support 
to  the  analysis  used  to  predict  the  linewidth  anticipated  from  the  proposed  high  power 
arrays. 

.Arrays  with  more  that  1  mW  of  power  would  contain  thousands  of  junctions  each 
requiring  a  bias  current  of  tens  of  milliamps.  For  the  arrays  described  above,  the 
junctions  were  biased  in  parallel  so  that,  for  an  N  junctions  array,  the  bias  current 
would  be  N  times  the  bias  current  of  a  single  junction.  This  arrangement  was  clearly 
impractical  for  very  large  arrays  although  it  had  the  advantage  of  permitting  a  greater 
scatter  in  the  junction  parameters  i.  e.  and  R.  .An  important  part  of  this  contact 
was.  therefore,  the  demonstration  that  successful  operation  could  be  achieved  in  arrays 
in  which  the  junctions  were  biased  in  senes  so  the  total  bias  current  was  equal  to  that  of 
a  single  junction.  We  designed  such  an  array,  based  on  the  quasi-lumped  design,  which 
was  fabricated  at  IBM  using  the  newly  developed  PARTS  technology.  Since  this  array 
was  a  small  part  of  a  much  larger  wafer,  the  values  of  J ^  and  R  were  dictated  by  the 
requirements  of  other  circuits  and  were  far  from  optimal  for  our  sources  so  that  we  were 
only  able  lo  design  for  a  power  level  of  1  pW.  The  resulting  100  junction  array, 
however,  worked  perfectly  delivering  the  1 /iW  t)f  power  to  a  5011  load  on-<hip  at 


300  GHz.  This  was  the  first  demonstration  that  the  series  biased  design  would  in  fact 
perform  as  predicted.  These  results  were  presented  at  the  1992  Applied 
Superconductivity  Conference  [5]  and  are  described  in  detail  in  appendix  E. 

The  final  phase  of  our  work  under  this  contract  was  to  attempt  to  combine  these 
various  developments  to  demonstrate  a  submillimeter  wave  source  with  100  /iW  of 
power.  This  500  junction  source  was  designed  using  the  distributed  array  concept 
demonstrated  above  (App.  C).  It  was  also  fabricated  at  IBM  using  the  PARTS  process- 
-in  fact  it  was  the  very  last  circuit  fabricated  before  the  line  was  shut  down  as  part  of 
the  IBM  retrenchment.  This  time,  howev  u,  we  were  able  to  include  the  source  as  part 
of  a  wafer  with  design  parameter  much  closer  to  the  ideal.  The  resulting  .source  again 
operated  as  predicted.  A  maximum  power  of  47  fj,W  was  delivered  to  a  68  H  load  on- 
chip  at  394 GHz  with  as  much  as  10 /rW  (still  higher  than  any  other  report  for  a 
Josephson  effect  source)  observed  at  500  GHz.  The  factor  of  two  difference  between  the 
designed  and  measured  power  is  accounted  for  since  the  actual  values  of  J ^  and  for 
the  fabricated  wafer  differed  somewhat  from  those  expected.  Using  the  measured 
junction  parame'iers.  the  observed  power  demonstrated  fully  coherent  operation  of  the 
source  so  we  declare  the  test  to  be  a  success.  V'arious  aspe  cts  of  this  work  were 
presented  at  the  1993  Terahertz  Technology  Conference  at  UCLA  [6]  and  at  the  1993 
ISEC  in  Boulder  [7]  (appendixes  F  and  G  respectively)  and  have  recently  been 
published  in  .Applied  Physics  Letters  [8].  This  work  is  presented  in  detail  in  .Vpp.  H. 

Two  Ph.  D.  theses  resulted  from  work  under  this  contract;  "High  Frequency 
Wave  Sources  Using  .Jo.sephson  .Junction  Arrays"  by  Kelin  Wan  in  1991  imd 
“‘Development  and  .Application  of  Josephson  Junctions  as  Submillimeter  Wave  Sources" 
by  Baokang  Bi  in  1993. 
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REFRACTORY  SUBMILLIMETER  JOSEPHSON  EFFECT  SOURCES 


K.  W«n,  B.  Bi,  A.K.  Jain,  LA.  Beaer,  S.  Han.  VV.H.  .MaUison,  and  J.E.  Lulcens 

Department  of  Physics 
State  University  of  New  York  at  Stony  Brook 
Stony  Brook.  New  York  11794 


ABSTRACT 

Niobium  Joeephson  effect  array  osdll&ton  have  been 
fabricated  and  tested.  These  derioes  deiiver  about  1  iiW  oi 
power,  in  the  inbmiUiffleter  wave  raage,  to  30-60  fl  load 
resiston.  The  preseat  upper  frequency  Limit  of  about  SOO  GHs 
appears  to  be  set  by  losses  id  the  supercouductiag  microstrip. 
Data  on  the  submilJiffletet  wave  su.'face  impedance  of  the 
Biobiuffl  microstnp  and  details  of  the  tnlayer  junction 
processiDg,  which  fires  value  of  at  2K  up  to  800  mV,  will 
be  presented. 


OfTRODPCnOW 

The  interest  in  submillimeter  wave  generation  unng 
Josephsoa  junctions  and  Joscphson  junction  amys  has 
inaeaied  in  the  pact  few  yean  because  of  the  possibility  ai 
using  them  as  local  osallators  for  SIS  mixen,  in 
communication,  and  radio-astronomicni  observation,  etc.,  das 
to  the  scataty  of  other  solid  state  sooree  ta  this  frequency 

ruga  A  Josephsoa  juction  is  a  voltaec  control  oscillaior 
(VCC)  with  the  potutisl  for  wide  tunabiLty  ud  a  fui  tuning 
rate.  Three  types  of  Josephsoa  effect  osoUaton  have  ben 

demonstrated  ,  namely,  the  resonut  Ouns  oscillator  ,  the 

I  C*4 

Ou  flow  osollator  ,  ud  small  junction  amys.  The  flitt 
two  lypet  nse  long  Josephsoa  junions  w.ih  the  juction 
dimennoa  much  latici  ihu  the  Josephsoa  puetrauoa  depth. 
In  this  paper,  we  will  focus  on  the  smail  iuncuoa  arrays. 

To  overcome  the  disadvutsgcr  of  a  stagle  metioe 
loaroe,  such  as  tow  impeducc,  low  output  power  icvel,  and 
large  liarndth,  u  array  having  a  large  number  (N)  of 

juctioni ,  all  of  which  are  phue-loched.  is  employed.  The 
output  power  irvsd  of  such  u  amy  lato  a  matUMl  load 
inereases  as  N  ,  while  the  Uuwidth  of  the  rsdisUoa  suited 
from  artsy  deorcsses  at  Isut  u  1/N.  Tb  trtu  the  amy  as  a 
lumped  arctui,  all  of  tee  juctioas  in  the  amy  ihonJd  be 
placH  within  about  one  eghth  of  a  wavdugih.  Thne  ,  tte 
number  of  juactious  is  liuied  br  the  unimum  spacing 
between  junctions  (shout  10  to  subsiuiially  less  thu  the 
optimum  number  of  junctiou.  To  arenmvat  tUs  problem, 

the  concept  of  s  diitnbuted  amy  tms  introduced.  Here  all 
junctions  arc  placed  one  waveJength  apart  so  that,  as  the 
operating  frequucy,  ihc  array  looks  like  a  lumped  orcmi. 

Previously  ,  we  have  demonstrated  that  disinbuted  amp 
of  Joscphson  junctions  having  40  junctions  could  delivar  1—7 
liW  of  power  onto  20-100  0  load  resistors  in  the  freqeeaey 
ruge  of  330  GBs  is  4S0  GHi.  The  upper  Umit  wu  set  by  the 
onset  of  the  large  locses  m  ihc  lesd  alloy  superoonduciing 
miaDstnp.  Those  arrays  were  made  of  Cu-thunted  lead  alloy 
tunnel  junctions  ud  suifeied  from  the  u^uaJ  problems  of 


Nfanusenpt  rrceved  September  24,  i09Q. 


non-uniformity  ud  cyclability  characteristie  of  this 
technology.  In  this  paper,  we  report  our  recut  remits  on 
similar  amys  with  An-shutad  Nb/AlOz/t^  juctions  , 
along  with  the  fabrication  details  ud  characteniues  of  t^  sU 
refraaory  Nb/AlOs/Nb  juctions.  Ideasuremuti  on  the 
siirfsce  impeduce  of  the  niobium  films  ,  unng  a  resonator 
drivu  by  a  resistively  shuted  Joscphson  juction  source  ,ate 
slio  diantssed. 


Nb/AlQi/Nb  TimNtL  JUNCTIONS 

1.  f  rmuhiiOB 

Juctions  were  fabrieatsd  u  outlined  below  osuii  the 
very  luccesstui  ntobinm  tnlayer  technoio^  (with  some  locnl 
venations)  which  has  been  developed  encniiveiy  dtutng  the 
last  de,'nde.‘''>4 

The  system  used  for  tht  Nb/AlOs/Nb  depocitios  in  the 
presut  work  coutitcd  of  multiple  sources  indndinr  s  DC 
magnetron  sputtenng  gu  for  Nb  depositiOB  ud  a  RF  inverted 
cylindficai  sputtenng  gu  tot  A1  depositioe.  The  miem  wu 
cryopumped.  having  a  base  preutut  below  2x10*'  Terr.  The 
ssmplas  were  mouted  on  a  water  cooled,  rotating  stage  which 
could  be  moved  over  the  vanou  sourcee.  The  distance 
betweu  the  sample  ud  the  target  wu  approismateiy  5  inchu. 
This  ruhu  Isrge  tpaoag  wu  used  .along  with  cooling  waut, 
to  euurt  the  substrate  temperuure  did  not  rue  dunni  the 
dcpoeiuoe  ud  also  to  reduce  the  stitu  in  the  ilms.  wc  feud 
that  prevnttng  substrate  heetieg  during  deposiura  wu  ouaal 
to  obiau  high  uality  juctiou. 

The  Nb/AiOi/Nb  txtlayer  tt'cnctazc  wu  putened  by 

muu  of  lift-off  uang  a  PVMA-AI-PMMA  tii-levcl  muk* 
rather  ibu  the  more  soavniional  method  of  etching  a 
ujuformJy  coatiid  tobiitatc.  The  lift-off  method  kept  the  film 
area  smsil,  rsduong  the  itieu  in  the  film.  This  wu  speoally 
importui  for  high  qnaiity  small  wea  juctiou.  The  tn-level 
muk,  iutcad  of  sutgi*  le^  math,  provided  e  better  tapered 
cdg^rofilc  thu  the  smelt  level  math. 

The  Nb/AJOxjNb  iruayer  itracturc  wu  made  by 
leqnuual  depoation  in  the  same  chamber.  The  1^  bau 

eltareds  wu  about  1000  X.  The  A1  vansd  from  10  A  to  400 
X  .  And  the  Nb  coutm  electrode  wu  aboei  800  X. 

Zvaporution  rates  tot  niobium  ud  sluminum  were  7-8  X/sec 

and  3.S-3  X/iee.  respectively.  The  tuuel  bsmer  wu  formed 
by  thermsJ  oadaiioe  of  the  Al.  Several  mesu  of  osodstma, 
iaduding  flowing  ozygu  gu  coniroUtd  by  mass  flow 
onatrollci.  otcygsn  pluma,  sad  10  %  of  Oi  mued  with  SO  %  of 
Ar>.  were  ineu.  It  wu  fe  und  ihsi  the  psntsl  pressure  of 
ozygu  deiennioet  the  osdc  thickacss  if  the  ondstioe  time  is 
kept  coutui.  To  vary  the  cxitical  currai  deasity,  tu 

pressure  of  the  ozygu  wu  dhed  u  5  mTotr  and  the  ondatioa 
time  wu  vanco  boot  several  miastat  up  to  one  hour.  This 


Mve  cnticai  cuneai  dusines  *sryiag,  u  shown  la  Fig  la, 
uoo  several  hundred  to  about  30  iiA/rm>.  where  the  gap 
voltage  of  the  junction  wu  reduced  to  about  2.5  mV 

(nu.'MA4«l«30O-3359$OI  (in  *5  iwi  IEEE 
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(»)  Jt  TS.  mbdttion  time.ioi,  the  mien  »how*  the 
Tiji»uon  of  \  ,  with  loi;  (b)  ]— V  di*r»aeiiit)c  of  k  high 
qnihtv  junction  with  V,w800  mV;  the  injeri  »how*  the 
J— V  chaxictenstic  of  jnnctioc  with  high  critic^  current 
dentity  showing  the  gap  roltige  abOTc  2.i  mV. 


The  pattern  of  the  junction  counter  deetrode  wa*  defined 
by  elearon  beam  lithography  using  a  tingle  layer  FMMA 

with  a  thickness  from  5000  X  to  1  pm  ,  depending  on  the 
of  the  junaioc  needed.  The  sample  was  etched  by  reactive 
ion  etching  (RIE)  with  CFi  gas  mixed  with  15  %  of  oxygen. 

1200  A  of  SiO  was  then  deposited  through  the  same  mmk  to 
form  an  insulating  layer,  isolating  the  counter  electrode. 
During  the  SiO  deposition,  the  samples  were  tilted  about  20 
degrees  and  rotated  to  that  the  edges  of  the  junction  counter 
electrode  could  be  fully  covered  and  a  uniform  SiO  layer 
formed  A  niobium  wiring  layer  was  deposited  at  the  aame 
rate  as  the  Nb  for  the  triiayet  Before  the  winne  depootion, 
the  surface  of  the  counter  elearode  was  cleaned  by  a  RT  Aij 
plasma  lor  5-7  minutes  to  remove  the  Nb  oxide. 


2.  CharaetenratiOB 

The  I— V  characteristics  of  the  Nb/AlOi/Nb  junctions  show 
that  the  junctions  can  be  of  high  quality  uiduding  high  V», 
high  gap  voltage,  and  a  narrow  gap  transiiion  width  Fig. lb 
shows  the  I— V  curve  of  a  junction  with  the  optinuxed 

labncation  parameter  using  15  Jl  AJ  and  oxidation  time  of  SO 
aunnies  This  gave  a  V,  of  about  800  mV  and  Jf  of  1500 
A/cm*  at  2  K  The  insert  of  Fig. lb  shows  the  I— V  curve  for  a 
junction  with  high  cntical  current  density  (40  kA/cm>) 
Although  the  junaion  has  large  lubgap  leakage,  its  gap  is  still 
above  2.5  m\'  Junctions  with  high  cntical  current  density  ire 
most  desirable  fo.'  Josephsou  efiect  oscillators,  while  the 
leakage  is  not  imponsni  since  the  junctions  are  to  be 

fCtalllVCjT  shufiltfd  )& 

Although  levera;  tens  of  A  of  A1  weie  geneisllv  used  to 
make  sure  the  niobium  surface  was  well  covered  by  Al,  it  had 


been  shown  by  X-ray  photoemission  ipectroecopy  ”  ,  and 
anodiaauon  spectroscopy  ,  that  only  pan  of  the  Al,  i.e., 

0—10  A,  is  oxidised  The  ecira  layer  of  unoxidixed  Al 

can  cause  a  reduction  in  the  junction  quality.'*  ”  To  examine 
this,  the  thickness  of  Al  was  vaned  ayitematicaliy  and  the 
dependence  of  the  junction  quabty,  V,,  on  the  ttuckness  of  Al 
vw  measured  (Fig.  21  The  junction  quality,  as  measured 
both  by  V  ,  and  the  aharpnets  of  the  gap,  is  a  snaromum  lor  15 

A  of  Al,  decreasing  steadily  for  both  greater  and  leaser 
thicknesses.  These  results  are  conautent  with  those  obtained 

by  Buggins  and  Gnrvitch. 
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F'sg.2.  Dependence  of  the  junction  quality  parameter  V» 
and  the  transition  width  of  the  gap,  AVmeasured  at  2K 
on  the  Al  thickness. 


For  Bie  in  Josephson  junction  array  oadUaton,  the 

junctions  are  required  to  have  a  small  fic  (Arw2irCIrEj/^t)  and 
an  1-V  charactrciitic  that  is  non  hysterestic.  This  was 

achieved  using  a  gold  shunt  resistor  with  a  30  to  50  A  Cr 
underlayer  for  adhesion,  M  shown  schematically  tn  Fig,3a 


(») 


I 

F^fi.3-  (i)  Schem*iic  of  the  AihHihuniH  Kb/AJOz/Nb 
j^ctioc;  (b)  J— curve  of  the  Aii-tauried  Nb  junction, 
the  in*en  u  the  equjviJeiil  orcuji  of  Ibe  junetjon, 
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Tjrpicailj,  300  X  of  Au  vfia  ’osea,  ji'nng  a  tesjiiaiiM  aboat  1  Ci 
pec  squixe.  Fig, 3b  aiso  ov«  tbe  i-'.'  caaractmsuci  of  the 
An-ihiuited  Nb/AlOi/No  junction  along  with  lU  eqaiTtleot 
dtcoit. 

The  junctionj  used  in  the  subnoiUimeter  w«Te  generator 
are  2x3  inn^  shouted  bj  an  approxiniaielx  2  jaq  long  leiiftor. 
The  critical  currents  of  the  junaions  vary  from  O.S  mA  to  1 
slA.  The  ST>ea£c  capaatance  ,  c  ,  of  the  nnctioni  wai 
determined  from  the  zero  field  step  (ZFSi,  giving  e  ;  40 
lF/;im>.  The  inductance  of  the  junction  shunt  resitor  loop  srae 
enimated  to  be  about  L  ^  0.1  pH.  These  parameters  are  tmail 
•Bongh  for  stable  operation  of  the  submiilimetor  sraTC 
generator  (i.e.,  A,  ^^2tLIc/^o<<1)- 


^n)/A10^/^fb  JUNCTION  ARRAYS 

The  design  of  Nb/AlOi/Nb  junction  array,  ahovn 
schematically  in  Fig.  4a.  it  similar  to  that  of  the  lead  alloy 
jonctian  arrays  ,  which  were  diicusted  in  detail  in  onr  prenons 

works.  Becanse  of  the  non— uniformity  of  the  jonctian 
asdcal  current  and  the  shunt  resistance  (  about  lOK)  ,  a 
parallel  biasmg  scheme  ,  which  compensates,  to  the  first  order, 
m  this  Tarnation,  was  still  used  All  of  the  junctions  are 
spaced  to  be  one  wavelength  apart  (about  350  can)  at  350  GHa. 
A  letutor  (about  20  fl)  is  placed  at  one  end  of  the  array  to 
nimnlat*  (lu  load.  The  r/ current  flowing  through  the  load 
resiftor  is  monitored  by  a  detector  Josephson  junction  placed 
at  one  end  of  the  load  resistor.  The  array  and  load  ate  r/ 
grounded  by  using  a  quarter-wavelength  miaostnp.  Direct 
grounding  (microstnp  line  shorted  to  the  ground  plane)  hat 
alto  been  tried  with  siimlar  resulti. 


The  critical  enrrent  and  the  shunted  resistance  of  the 
junctions  used  for  these  arrays  were  0.5-1  mA  and  0.5  to  0.7 
n,  respectively,  giving  the  leRi  product  of  about  500  jiV,  and 
the  I-V  curve  of  the  array  is  non  hytteresuc.  The 
•trays  aeng  Nb/AlOi/Nb  junctions  generally  operate  over  a 
ranch  wider  frequency  range,  from  100  GHs  to  above  500  GHi, 
and  exhibited  increased  cydibUity  and  parameter  uniformity. 
The  arrays  can  deliver  a  power  of  about  1  nW  to  20-60  fl 
load  resistors  u  seen  in  Fig.  4b,  showing  the  power  to  the  load 
vs.  frequency  u  the  amy  bias  is  varied  The  power  level  was 
determined  by  using  the  Josephson  junction  detector  by 

measuring  the  amplitude  of  the  first  Shapiro  step.'**  The 
insert  of  Fig.  4b  shows  the  details  of  this  power  variation  n>ar 
280  GHs.  One  can  see  that  the  power  vanes  with  a  penodiaty 
at  about  4.5  GHs  which  corresponds  to  c/2U,  where  e  and  U  are 
the  phase  veloaty  and  the  total  length  of  the  array.  As  can  be 
seen,  appreaable  power  ii  obtained  from  the  array,  not  only  at 
demgn^  frequency  v,  and  ,  (where  the  junctions  are 
spac^  by  A  and  A/2,  respectively)  bnc  also  over  a  wide  range 

of  freetnencies  satisfying  fo  •  n  c/2k  (n  an  integer).  A  theory 
far  this  breed  band  operation  has  yet  to  be  worked  ont. 

We  have  not  yet  made  a  direct  measurement  of  the 
Unewidth  of  the  radiauon  from  this  array.  The  RSJ  analysis, 
which  agrees  within  a  factor  of  two  with  the  measured 

Unewidth  far  paralld  arrays  in  the  mioowave  range  was 
used  to  estimate  the  Unewidth,  far  these  snbnullimetet  arrays. 
The  Unewidth  of  the  radiation  at  i/g  is  thus  estimated  to  be  60 
KHs  at  2  K.  and  130  KBi  at  4.2  K,  respecPveiy  from  the 
measured  differential  retiitance  ba^  on  this  lumped  oremt 
inalysu. 


FIg.5.  Schematic  of  tlie  resonator  driven  by  a  Josephnn 
junction  sonree  u  used  for  surface  impedance 
meunrementa. 


SERFAa:  FiESISTANCE  OF  THE  NnCHOSTRIP 
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ngA  (a)  Schematic  of  the  parallel  biased  distnbuted 
array  The  junctions  (xi  are  placed  ai  350  urn  iniervais 
along  the  serpentine  micipsinp  The  load  lesistor  and 
detecior  junction  are  shown  in  ihe  upper  leli  The  daihed 
Unes  show  the  erte.ni  ol  ihe  ground  plane  under  the 
microstnp.  (b)  Power  vs  fresiuencv,  the  insert  is  ihe 
expansion  near  260  C ITi ,  ho wmg  a  pfrioiimi  v  of  lOont 
4.5  GHi. 


In  order  id  measure  the  phase  veloaty  and  lots  of  th«i 
iupertondncimif,  microiinp  in  the  enbmiiUmeter  wave  range,  a 
microstnp  resonator  which  uses  a  Josepheon  junction  u  both 
the  source  an/i  detector  hu  been  deveioped.  The  tchemaaic  of 
this  resonator  is  shown  in  Fig. 5.  Our  approach  differs 

somewhat  from  previous  measurements  of  this  type  in  that 
a  motuvely  ihunted  junction  has  been  used  in  the  resonator. 

In  the  absence  of  r/ currents  flowir,g  in  the  resonator,  the  1-V 
curve  of  the  juction  is  given  analytically  by  the  weil  knows 
RSJ  model.  The  biu  lea^,  which  are  placed  next  to  the 
junction  aa  as  a  high  imp^ance  transmisBon  line,  which  tf 
properly  terminate^  limply  causes  a  tUght  shift  in  the  thnnt 
reDitance. 

The  rfcvmnu  which  flow  in  the  microitnp  when  the 
wavelength  A  e  2  4  where  /  it  the  length  of  the  resonator,  aa 
u  a  perturbation  on  the  RSJ  loloiion  The  junction's 
charaaeniiii  frequency,  given  by  its  URj  prodna,  is  selected 
to  be  near  or  below  the  fiequenaa  of  Intcreat.  For  these 
conditions,  the  well  estsblithed  penurhauon  spp  'ch  agrees 
very  weU  with  direa  computa  timulaiions.  The  j.maion'1 
response  it  thus  weU  chatacienied  and  euy  to  aiialy-. e.  The 
mioDg  of  Josephson  osaUstions  in  the  junction  with  .  he  rf 
micioiiiip  currents  prwlncei  s  ihifl  in  the  biu  current 
required  for  a  given  opealing  freouency  or  volisge.  Thii  ihift 
1/ II  given  by  if  «  const  Re.(  1/Z)  where  Z  it  ii.  r/  loop 
impedance  of  the  itrorture  ihown  :o  Fig  5  Thus  ihe  difTcience 


between  the  jonctioc'c  bias  cnrresii  with  and  witbotii  tbe 
nurroiinp  resoDance  ic  just  a  reconast  peak  lu  nbown  id  Fig.6 
Here  the  data  with  the  best  fit  BSJ  tuckfnoacl  gnbttacted  are 
ahowD  as  points,  while  the  best  fit  to  the  theo77  i<  the  solid 
lute. The  resonant  ireqnenc;  (IM  GEx)  prei  the  phaie  eelodt}' 


V(mV) 


ng.6.  Detail  of  letonance  peak  S/(m)  with  the  best 
theoretical  fit  (— ),  the  tnaert  is  the  I— V  ecrrc  of  the 
jonmon  in  the  resanatoi. 


V^»0.4r^  and  the  Q  (;  SO)  prm  the  lou  in  the  micrQStxip.  For 

the  micTOctnp  dunensioni  of  SOO  imt  loot  bj  10  ion  vide,  this 
Q  uspbes  a  ssriace  resstanoe  of  3  ntf)  for  the  lute  at  iJlK 


C0^’CLPSlg^ 


We  hare  diicossed  the  (abricaticis  and  characierisaticin 
cf  Nb/AlOz/Nb  junctions  haYui{  Va  of  up  to  lOOmV  at  2K 
Arrant  osui£  Nb/AlOz/Nb  junctions  hare  been  (abnealed  and 
tested  These  am.rt,  which  deiieeted  about  1  pW  of  power  to 
30-60  fl  load  recitors  operated  over  a  wide  fi^uenc;  r^e, 
from  100  GEt  to  SOO  G^.  and  exhibited  improved  c^buitj 
and  parameter  Ttnilormitr  oompared  to  arra.n  nsin£  lead  411(7 
junctions  The  surface  impedance  in  the  lubmiUimeior  wave 
rbsge  was  measured  ustn;  a  misostnp  le^naioi  driven  I7  a 
Josephson  junction,  giving  c  3  mfl  at 
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Josephson-Junction  Driven  Submillimeter 
Wave  Microstrip  Resonator 

Baokang  Bi,  K.  Wan,  Wenxing  Zhang,  Siyuan  Han,  and  James  E.  Lukens 


Abstract  — A  microstrip  resonator  driven  b)  a  resistively 
shunted  Josephson  junction  has  been  used  for  submillimeter 
wave  surface  resistance  measurements  to  about  400  GHz.  A 
simple,  analytical  analysis  of  Ibe  resonator  induced  structure  in 
the  Junction's  /-f  characteristic  gives  the  resonator  surface 
resistance  in  the  range  of  2-160  mO  with  an  accuracy  of  better 
than 

1.  Introdlction 

HE  DEVELOPMENT  of  thin  film  devices  in  the  submil¬ 
limeter  wave  range  has  become  an  active  research  area, 
leading  to  ihe  need  to  understand  the  physical  properties, 
such  as  the  surface  resistance  and  of  materials  in  this  fre¬ 
quency  range.  Below  about  200  GHz,  resonant  cavities  cou¬ 
pled  to  solid  state  (e  g  ,  Gunn)  oscillators  are  useful  for  these 
measurements,  while  for  /  a  1  THz.  infrared  techniques 
become  available  For  intermediate  frequencies,  due  to  the 
lack  of  suitable  sources,  surface  resistance  measurements 
have  proven  more  difficult  One  choice  for  these  frequencies 
is  to  use  superconducting  Josephson  junction  oscillators.  The 
Josephson  junction  is  a  natural  voltage  controlled  oscillator, 
with  an  oscillation  frequency  related  to  its  dc  bias  voltage  by 
/=  2el'  /i  =  48.^  GHz  mV  For  typical  junctions  (e  g  , 
Nb 'AlO^  Nb)  oscillations  are  sustained  to  over  1  THz. 
When  coupled  to  a  resonant  structure,  the  Josephson  junction 
radiates  into  the  resonator,  while  the  resonator  current  inter¬ 
acts  with  the  junction  The  results  of  this  interaction  are 
reflected  in  the  junction  dc  /- !■'  curve,  which  carries  the 
information  about  the  resonator  This  phenomenon  has  been 
studied  theoretically  by,  eg,  Likharev  (1) 

There  are  two  useful  systems,  distinguished  by  the  mode 
of  coupling  to  the  resonator  The  first  one,  which  has  been 
widely  used,  consists  of  a  weakly  damped  SIS  junction 
coupled  to  a  parallel  type  of  resonator  This  system  has  been 
used  to  improve  the  impedance  matching  of  SIS  mixers  [2] 
and  to  measure  the  phase  velocities  of  dielecmc  materials 
!3),  (4)  The  second  system  is  made  of  a  highly  damped 
Josephson  junction  coupled  to  a  senes  type  of  resonator  The 
resonances  in  this  latter  system  have  been  experimentally 
demonstrated  by  using  point  contact  (3)  and  nncrohridge  (6) 
junctions 

in  this  paper,  we  report  work  on  deseUiping  a  resistively 
shunied  Josephson  junction  (RSJ)  dnsen  subniillirneier  wave 
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microsirip  resonator,  as  shown  in  Fig.  ).  This  enables  us  to 
measure  both  the  surface  resistance  /?,,  defined  as  the  resis¬ 
tance  per  square,  and  the  relative  dielectric  constant  e,  in  this 
frequency  range.  The  interaction  of  the  RSJ  with  the  res¬ 
onator  shows  up  as  a  bump  (see  Fig.  2(c))  on  the  smooth 
/-  V  curve  of  the  unperturbed  junction  and  can  easily  be 
calculated  by  using  either  numerical  analysis  or  perturbation 
analysis  [7]  starting  from  the  analytical  RSJ  I-V  curve.  The 
resonance  due  to  the  coupling  between  the  RSJ  and  the 
re.sonator  has  two  main  features.  The  first  feature  is  the 
center  resonant  frequency ,  which  depends  on  the  dielectric 
material  used.  The  second  one  is  the  resonant  width  or  Q 
factor,  which  is  related  to  the  loss  in  the  system.  By  using  the 
RSJ  as  both  the  submillimeter  wave  source  and  the  detector, 
both  of  these  two  features  can  be  readily  obtained  from  the 
junctions’s  dc  /-  ^  curve.  Therefore,  the  values  of  /?,  and 
t,  can  be  extracted.  Other  advantages  of  using  the  RSJ,  such 
as  a  smooth  /-  V'  background  curve  and  continuously  tunable 
operating  frequency  are  also  important  considerations  for  its 
choice, 

II  Device  Prepailation 

Fig  1  shows  the  structure  of  the  microstrip  resonator.  A 
patterned  ground  plane  is  deposited  first  The  supponing  SiO 
dielectric  is  then  thermally  evaporated  in  a  chamber  with 
base  pressure  in  the  10  ’  lorr  range.  Finally,  the  microstrip 
transmission  line  with  the  RSJ  at  its  center  is  fabricated.  The 
microstrip  transmission  line,  the  SiO  dielectric  layer,  and  the 
ground  plane  form  an  open  ended  microstnp  resonator.  The 
materials  used  for  making,  the  resonator  are  indicated  by  the 
noution  MS  SiO  GP  Typically ,  the  thickness  of  SiO  dielec¬ 
tric  d  is  700  nm,  and  the  width  of  microstnp  w  is  10  fjtn 
The  length  of  resonator  /  is  determined  by  the  desired 
frequency  and  the  phase  velocity  in  the  dielectric,  for  exam¬ 
ple,  /  =  200  f^m  for  a  3(X)-GHz  fundamental  resonance  with 
thennally  evaporated  SiO  dielectric  The  RSJ's  used  here  arc 
mainly  Au  shunted  Nb,  AlO,  /Nb  junctions,  using 
STv/AlO,  NT)  trilayer  technology,  as  given  in  [8)  Pb-alloy 
junctions  are  also  used  occasionally  and  are  descnlsed  in  [9] 

One  important  design  parameter  of  the  RSJ  is  its 
product,  where  /,  is  the  critical  curreni  and  is  the  shunt 
resistance  The  spectrum  of  Josephson  oscillations  is  rich  in 
higher  order  harmonics  [7)  when  the  bias  voltage  is  low 
()'</,,/<^).  To  ensure  mon(Khr,imaiic  Josephson  oscilla- 
non,  should  be  smaller  then  the  designed  rcsonani 

frequency,  /,  <  /j/q  '2f,  On  the  other  h.ind,  the  junction 

resistance  should  be  kept  small  for  the  purposes  of  both 
reducing  the  resonator  loss  and  suppressing  the  subgap  struc' 
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•  ,tl!  'it* 

Fi|  1  Siructure  of  rrucrosinp  reson^ior  RSJ  is  iocatfii  4i  the  center  of 
rT\jcro5tnp  /-  y  curve  u  uken  b)  using  4-iemrunaJ  meaiuremenis 


Fig  7.  («)  Ei^uivaJeni  elcriricaJ  circuH  for  a  junction  coupled  to  a  res¬ 

onator  ft-L-C  senes  resoneni  circuM  is  used  to  simutaic  the  microsinp 
resonator  for  iheorciicil  calculation  fb)  High  frequenc)  equivalent  circuit  of 
(a)  RSJ  IS  nrodelcd  as  a  RF  oscillator  »iib  internal  source  resisuiKc  Rj. 
The  OKillaior  produces  a  RF  current  in  loop  (c)  Comparison  of 
calculated  results  from  numencat  anaJvsis  (solid  lirse)  and  perturbation 
theory  (solid  square)  Unperutrbed  /•V'  curse  (dash  line)  is'deftned  as 
background  /-  L  curve  The  bump  in  the  curve  is  due  to  the  interaction 
vkith  the  resonator 

ture,  which  is  essential  to  obtaining  a  smooth  background 
I-  y  curve.  The  ground  plane  is  only  about  10  *im  wider 
than  the  microstrip,  so  that  the  bias  leads,  placed  close  to  the 
RSJ.  can  come  off  the  ground  plane  immediately.  This  re¬ 
duces  the  RF  leakage  through  the  dc  bias  leads  to  a  mini- 
mum.  Table  I  summarizes  the  important  dimensions  and 
fabrication  conditions  of  the  materials  used  for  making  the 
microsirip  resonators. 

lU.  Analysis 

The  microsirip  re.sonaior  is  characterized  by  the  attenua¬ 
tion  consum  a,  the  wave  number  0,  and  the  characteristic 
impedance  Zq  which  is  determined  by  unit  length  inductance 
L  and  capacitance  C  of  the  microsirip  transmission  line, 
Zq  =  f/L/C)' TTie  resonator  has  length  /  and  width  w. 
The  tJiickness  of  the  dielectric  is  d  Fig  2(a)  shows  the 
equivalent  electrical  circuit  for  a  RSJ  coupled  to  a  resoi  ator. 


TABLE  I 

Fajiucatx)»*  CoNorrioNs  *nd  Physical  pAtAMrrziis  of  Matiiuau  Tutid 


Mat^ruJs 

PrepAraiion 

Thicknesj  (nm) 

»(0'  'm  ■ 4  2K 

Cu 

e-beani 

150  - 

3  9  «  to' 

Au 

C'bcam 

too  ■  ■- 

8  8  «  10’ 

Nb 

DC  spunered 

1  (THon  Ar 

300  -  2  6 

2  8  to  111  lOKi 

*45  mcjLsured  b)  using  SIS  runnel  junctions 


For  simplicity  in  numerical  simulation,  the  microstrip  res¬ 
onator  is  replaced  with  a  -C,  series  resonant  circuit, 

where  R,  =  o/Zq  represents  the  total  loss  in  the  microstnp 
resonator,  and  L,  =  IL/2  and  C,  =  2IC/r^  are  the  equiva¬ 
lent  inductance  and  capacitance  for  the  resonator  shown  in 
Fig.  1  at  the  fundamental  resonance.  The  effect  of  junction 
shunt  capacitance  is  neglected  for  our  sa.oples.  since  0^  m 
2el^R]Cjlh  <0.02.  Fig.  2(b)  shows  the  high  frequency 
equivalent  circuit  [7),  where  the  junction  is  mexJeled  by  a  RF 
oscillator  with  internal  source  resistance  Rj.  The  behavior  of 
this  RSJ,  when  coupled  to  a  resonator,  can  be  analyzed  in 
two  ways.  The  fundamental  method  is  by  the  numerical 
solution  of  the  set  of  differential  equations  describing  the 
circuit.  In  normalized  units,  /  =  ///f.  v  =  V jl^Rj,  r  = 
2el^Rjt/h,  these  can  be  wrinen,  for  the  electrical  circuit 
shown  in  Fig.  2(a),  as 

'6  =  'l(0  +  '2(0  (1) 

i|(0  =  '»(’■)  +  sins^Cr)  (2) 

=  IjCr)/-,  -i-  /S„/j(r)  0;,'/2(t)  (3) 

where  i^,  is  the  dc  bias  current,  »,(t)  is  the  current  flowing 
through  the  junction,  Ij(t)  is  the  current  in  the  resonator 
branch,  and  ifiir)  =  u(t)  is  the  voltage  across  the  junction 
r,  =  R,/Rj  is  the  nonnaJLzed  resistance  of  the  resonant 
circuit.  The  capacitance  C,  and  the  inductance  L,  of  reso¬ 
nant  circuit  are  included  in  the  parameters  0„  = 
2e/^RjC,lh  and  0,,  =  2?/,^,/^.  Equation  (2)  is  the  dy¬ 
namic  equation  for  a  resistively  shunted  junction  [1],  and  (3) 
is  the  differential  equa'ion  for  the  resonant  circuit.  Equations 
(l)-(3)  are  coupled  equations  and,  in  general,  can  only  be 
solved  numerically.  By  taking  the  time  average  of  ufr),  the 
dc  /-  y  curve  of  the  junction  is  obtained.  An  aJiemaiive  and 
more  instructive  method  for  understanding  the  b-tiavior  of 
the  RSJ  is  perturbation  theory  (7J.  In  this  theory,  the  interac¬ 
tions  between  the  RSJ  and  external  radiation  are  considered 
as  a  perturbation  to  the  RSJ  I-  V  curve,  which  is  given 
analytically  by  the  well-known  RSJ  model  (1).  'VMien  biased 
near  the  resonance,  the  junction  oscillation  produces  a  RF 
current  in  the  resonator  (see  Fig.  2(b)).  This  RF  current 
in  turn  mixes  with  the  junction  oscillation.  The  results  of  this 
interaction,  following  from  the  perturbation  theory,  are  re¬ 
flected  as  an  effective  dc  offset  i„  in  the  bias  current 

‘m{^)  =  ('•) 

where  ~  y„/l,Rj  is  the  normalized  magnitude  of 
Josephson  oscillation,  y  is  the  down  conversion  coefficient  of 
mixing,  and  yfu)  =  Rj  Re(l/Z)  is  the  real  part  of  normal¬ 
ized  admittance  for  the  loop  shown  in  Fig.  2(b).  The  offset 
current  has  been  normalized  to  /,.  Within  the  resonant 
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width,  both  >  and  can  be  treated  as  constants.  When  the 
junction  frequency  is  swept  across  die  resonance.  i„  -  u  has 
a  universal  Lorentzian  resonant  shape  with  its  Q  determined 
only  by  the  loop  admittance.  The  dc  /-  V  curve  is,  therefore: 

/(c)  -  /"(u)  -  /„(v).  (5) 

The  first  term  on  the  right  side  /“(c)  -  (1  +  c')'  '  is  the 
unperturbed  RSJ  I~V  curve,  which  we  call  the  background 
I~  V  curve.  The  second  term  on  the  right  side  is  the  perturba¬ 
tive  term,  !-V  curves  for  identical  parameters  have  been 
calculated  by  using  both  (l)-(3)  and  (5).  As  can  be  seen  in 
Fig.  2(c),  the  results  from  these  two  approaches  agree  well 
with  each  other.  Perturbation  theory,  however,  provides 
simple  analytical  results  for  the  structure  induced  on  the  I-  V 
cuive  by  the  resonator,  simplifying  the  dau  analysis 

For  the  resonator  used,  an  open  ended  miciostrip  transmis¬ 
sion  line  with  a  junction  placed  at  the  center,  tlie  loop 
impedance  seen  by  the  oscillator  is  the  sum  of  source 
impedance  and  the  transmission  line  impedance: 

Z  =  -f  2Zo  coth  (a -(- //3)|  (6) 

where  Zq.  the  characteristic  impedance  of  microstrip  line, 
can  be  calculated  using  the  equation  given  by  [10],  At 
resonance,  the  imaginary  pan  of  (6)  vanishes,  which  gives 
resonant  frequencies 


(>"  -  1)  t’p(<cft) 
2  / 


n=1.2.--  (7) 


where  /  has  been  taken  as  the  physical  length  of  the  resonator 
since  the  corrections  due  to  the  open  ends  are  calculated  to  be 
0. 1  )  is  the  phase  velocity  ,  the  formulation  of  which 

is  given  in  (10)  for  superconducting  microstrip  lines. 
is  the  effective  dielectric  constant  and  is  related  to  the  rela¬ 
tive  dielectric  constant  t,  by  =  (l/2)(«,  -f  1)  -(- 
(l/2)({,  -  1)(1  -t-  12df/H')“ '  ‘[1 1).  The  effect  of  microstrip 
dispersion  on  «,  is  negligible  here.  Assuming  Uiat  there  are 
no  other  losses,  the  surface  resistances  per  square  are  directly 
related  to  the  attenuation  constant  a  [11]  by 

f?.,  =  2awZo  (8) 


where  f?,,.  /?,;  are  the  surface  resistances  of  the  rnicrostrip 
and  the  ground  plane,  respectively  The  Q  of  resonance  can 
be  calculated  as 


Q 


1  ^0 

2  Rj  +  oIZq 


2  Rj  +  R,(l/w)  ■ 


(9) 


The  last  equal  .sign  assumes  that  R^  =  =  Rj;  if  the 

microstrip  and  the  ground  plane  are  made  of  the  same 
materials. 

For  data  analy  sis,  penurbation  theory  is  used  The  reso¬ 
nant  peak  is  obtained  by  subtracting  the  best  fitted  back¬ 
ground  i“(i')  from  the  expTCrimenial  /-  V  curve.  The  attenua¬ 
tion  a  is  then  extracted  by  fitting  the  resonant  peak  to  (4), 


where  the  loop  admittance  yti  )  is  calculated  from  (6)  and 
7i'„  is  simply  treated  as  a  scaling  constant  The  surface 
resistance  follows  from  (8l,  assuming  that  there  are  no  other 
sources  of  loss  in  the  microstrip  . 

IV.  Results 

A  typical  experimental  curve  from  a  NT),  SiO  Nb  resonator 
is  shown  in  Fig.  3  with  the  junction  having  l^R^  =  15'*  ,  V, 
Rj  =  96  mf),  and  tiie  nicrostrip  having  /,  w  =  30  and 
Zq  =  10  n.  The  strucnire  appearing  on  the  l-V  curve  in 
Fig.  3(a)  is  due  to  the  resonance  from  the  microsttip  res¬ 
onator.  Fig.  3(b)  is  the  resonant  peak  after  the  background 
/“(i>)  has  been  removed  The  best  fit  to  (4)  gives  Q  ~  54, 
/o  =  198  GHz,  implying  /?,  =  6.1  mfl.  The  effect  of  iher- 
trial  noise  is  negligible  for  this  resonance.  A  number  of 
samples,  made  of  different  materials  and  with  different  de¬ 
signed  resonant  frequencies-,  have  been  fabricated  and  tested 
Tfie  surface  resistance  values  obtained  from  those  expieri- 
ments  at  4.2K  are  summarized  in  Fig  4(a).  For  Nb  films,  the 
data  were  collected  from  several  different  resonators  made 
with  different  designed  frequencies  and  most  likely  having  a 
variation  in  R^  among  them.  The  solid  bne  is  the  calculated 
surface  resistance  of  the  NT)  rnicrostrip  transmission  line  by 
using  the  Mattis- Bardeen  theory  [12].  (13)  with  a  film 
thickness  of  /  =  300  nm.  a(lOK)  =  2  8  •  lO'D'  'm’ and 
2S  =  2.6  mV.  Presently,  the  experimental  data  cover  the 
frequency  range  up  to  400  GHz.  Phase  velocities  measured 
from  NT)  SiO/NTi  rnicrostrip  resonators  at  4,2  K  are  plotted 
in  Fig,  4fb).  Within  this  frequency  range,  is  almost 
constant.  By  using  the  average  phase  velocity,  taking  = 
83  nm  [14]  and  the  dielectric  thickness  d  =  700  nm,  the 
relative  dielectric  constant  for  our  thermally  evaporated  SiO 
is  calculated  to  be  t,  =  5.55  ±  0.14  which  is  in  good  ag- 
greement  with  data  published  for  frequencies  under  100  GHz 
[15],  The  solid  line  in  Fig.  4(b)  is  the  phase  velocity  calcu¬ 
lated  by  taking  «,  =  5.55,  where  the  complex  conductance 
from  the  Mattis -Bardeen  theory  is  used. 

V.  Discussions 

The  surface  resistance  values  obtained  from  die  previous 
experiments  have  been  attributed  to  conductor  losses  of  metal 
films  forming  the  rnicrostrip  resonators  However,  other  loss 
mechanisms,  such  as  radiation  loss  and  dielectric  loss,  also 
exist  in  the  resonator.  Accurate  calculation  of  those  los.ses  in 
our  resonator  is  very  difficult,  so  only  estimations  will  be 
given.  The  amount  of  energy  radiated  from  the  op>en  enos 
depends  on  the  design  parameters  and  operating  frequency  of 
the  rnicrostrip  resonator.  The  radiation  Q  factor  for  a  half 
wavelength  resonator,  such  as  ours,  has  been  calculated  by 
Lewin  [16],  For  Zq  =  10  fl,  /  =  200  GHz,  d  =  700  nm. 
and  f,  =  5.55,  one  has  Q,  =  13  000.  The  loss  tangent  for 
SiO  in  the  submillimeter  wave  range  at  low  temperature  is 
not  known  at  present.  From  the  available  data  at  4  2  K  and 
30  GHz  [17],  a  value  of  tandj  -  10“’  should  be  a  reason¬ 
able  estimate  for  tani^  at  several  hundred  gigahertz.  This 
gives  =  l/tan6^  ^  1000.  By  using  (9),  the  losses  from 
radiation  and  the  dielectric  can  be  expressed  in  terms  of 
surface  resistance.  Taking  //w  =  30,  then  R,,^  =  0.04  mfi, 
R^  =  0.5  mO.  Among  R,^  and  R^,  the  dielectric  loss  R^ 
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Fig  3.  (a)  Typical  expenmenial  l-V  curse  shovong  resonance  (b)  The 

resonant  peaii  (solid  square)  from  (ai  after  background  is  removed  Solid  line 
IS  the  best  fit  lo  penurbalion  theory  giving  Q  -  iA 


f  (GHz) 
(b) 


Fig  4  (a)  Esi'acled  surface  resistance  values  of  Au  (diamond)  from 

Pb  SiO  Au(l(X)  nm).  Cu  (iriangle)  from  Pb  SiO  Cu(150  nm).  and  Nb 
(circle)  from  Nb  SiO  Nb(300  nm)  ai  4  2K  from  SO  to  AOO  CHr  Solid  line  is 
the  surface  resistance  of  Nb  microstrip  ai  4  2K  vvith  2d  =  2  6  mV,  calcu 
laied  by  using  Maiiis-Bardeen  theory,  (b)  Phase  veiociiies  (square)  mea¬ 
sured  from  Nb  SiO  Nb  microstnp  resonators  at  4  2K  from  50  to  400  GHz 
Solid  line  is  the  phase  velocity  given  by  Matlts- Bardeen  calculation  with 
e,  =  5  55.  Noie  smal'  curvature. 


IS  ihe  largest,  and  could  be  as  large  as  0  .5  mO  With  those  in 
mind,  It  IS  sensible  to  conclude  that  the  dominant  loss  in  Fig 
4(a)  comes  from  the  conductor 

The  smallest  surface  resistance  thal  can  be  measured  b> 
using  this  methcxJ  is  limited  by  the  loss  frequency  voltage 
noise  and  by  the  total  internal  losses  including  oscillator 
shunt  resistance,  dielectric  loss,  and  radiation  loss.  The  total 
internal  loss  should  be  kept  as  small  as  possible  By  using  a 
suspended  ground  plane,  dielectnc  loss  could  be  reduced  to 
minimum  A  smaller  shunt  resistance  Rj  <  0  5A?,  should 
also  be  chosen  so  thal  the  loss  due  to  the  surface  resistance 
alstays  dominates,  as  shossn  in  Fig  3, 

The  upper  bound  of  the  measurable  surface  resistance  is 
primanly  determined  by  hosv  well  a  broad  resonant  peak  can 
be  recovered  from  the  background  I-V.  For  the  resonators 
made  of  superconducting  materials,  when  the  operating  fre¬ 
quency  />  2A  h,  the  surface  resistance  increases  dramati¬ 
cally  and  the  resonance  is  not  detectable,  this  limits  the 
highest  usable  frequency,  to  be  about  650  GHz,  if  Nb  is  used 
for  the  microstrip. 

The  remaining  uncertainties  of  the  suiface  resistance  data 
can  be  evaluated  from  known  parameters.  The  uncertainties 
of  width  Aw/w,  and  characteristic  impedance  Zq  are  about 
5  and  13%,  resjyectively,  as  estimated  from  fabrication  con¬ 
trol  over  the  microstrip  parameters.  Errors  of  Aa/a  from 
data  fitting  are  no  more  than  5%  .n  most  cases,  giving  an 
overall  uncertainty  of  deduced  surface  resistance  values  of 
about  25%, 

A  final  remark  on  our  experimental  results  is  that  the 
surface  resistance  values  are  measured  from  real  devices, 
which  have  gone  through  a  series  of  fabrication  steps.  Metal 
thin  films  niay  have  fine  grains,  stress,  and  surface  layers, 
especially  for  Nb  fil.ns.  All  those  factors  tend  to  deteriorate 
the  RF  properties  of  thin  films.  Therefore,  the  microstrip 
resonator  described  here  is  not  an  ideal  system  for  measuring 
the  intrinsic  losses  of  conducting  thin  films  unless  precautions 
are  taken  in  preparing  and  characterizing  the  thin  films,  but  is 
suitable  for  measuring  in  situ  physical  parameters  which  are 
very  important  for  designing  practical  %F  devices  I'kc  phase- 
locked  submillimeter  wave  oscillators  [7] -[9]  and  microstnp 
transmission  lines.  Placing  the  films  to  be  measured  on  a 
suspended  ground  plane  would  also  alleviate  the  problems. 
Work  on  this  i  chnique  is  in  progress. 

VI.  Conclusions 

In  summary,  RSJ  driven  submillimeter  wave  microstrip 
resonators  have  been  developed.  The  behavior  of  a  RSJ 
coupled  to  the  microstrip  resonator  has  been  well  character¬ 
ized  theoretically.  The  microstrip  resonators  have  been  used 
to  measure  the  submillimeter  surface  resistance  with  an  accu¬ 
racy  of  25%  as  well  on  the  dielectric  constant  up  to  400 
GHz.  This  method  can  be  readily  extended  to  a  higher 
frequency  range  by  using  superconductors  with  a  larger  gap 
energy  2A. 
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DISTRIBUTED  JOSEPHSON  JUNCTION  ARRAYS  AS  LOCAL  OSCIl,LATORS 

Banking  Hi,  S.  Han,  J  E,  Lukftis.  and  K.  Wan* 

Drparinifnt  of  Physirs 
at  Stony  Brook.  Stony  Brook.  N'\'  11794 


Atifiraci  -  The  power  and  linewidfh  of  the  radiation, 
near  250  GHz,  of  sriiail  distributed  Jorephson  junction 
arrav  sources  have  been  measured  using  a  single  chip 
integrated  source  and  receiver.  The  one-dimensional 
arra>s  were  of  a  new  design  w-ith  junctions  placed  in  !/1 
travfleog'h  lumps  separated  by  one  half  wavelength. 
The  meaisured  power  (  rr2p\V)  and  linewidth 
(~  10  MHz)  for  the  10  junction  sources  were  in 
reasonable  agreement  with  those  obtained  from  computer 
simulations. 

I.  INTRODUCTION 

Over  the  last  decade.  Josephson  junction  arrays  have 
been  demonstrated  to  be  useful  for  generating  radiation 
in  the  millimeter  and  subrnillimeter  wave  ranges[l-4]. 
The  first  generation  of  Josephson  junction  arrays  were 
lumped  circuits,  requiring  that  the  distribution  of 
junctions  should  be  within  one-eighth  of  a  wave  length 
a*,  the  operating  frequency.  Such  a  requirement  limits 
the  maximum  number  of  junctions  that  can  be  accom¬ 
modated,  and  therefore  the  maximum  available  jxiwer 
level.  The  “quasi-lumf>ed"[4]  array,  in  which  junctions 
were  placed  one  wave  length  apart,  was  developed  to 
remove  this  limitation.  Expierimental  results  on  quasi- 
lumped  arrays  have  demonstrated  their  ability  to  deliver 
higher  output  power  to  a  high  impedance  load. 
However,  the  design  of  the  quasi-lumped  array  results  in 
»  large  structure,  which  is  less  compact  and,  more 
importantly,  has  significant  transmission  line  loss  at  high 
frequency. 

In  this  paper,  we  report  an  improved  design  for  one 
dimensional  Josephson  junction  arrays,  incorporating  the 
concepts  of  both  lumpied  and  quasi-lumped  arrays.  In 
this  new  design,  junctions  are  separated  into  groups  with 
sn  inter-group  spacing  of  half  a  wave  length,  see  Fig.  1. 
This  structure  approximates  that  of  the  quasi-lumped 
array.  However  a  group  of  junctions,  which,  in  principle, 
can  be  distributed  over  a  length  up  to  1/2  A,  replaces 
each  single  junction  used  in  the  previous  design.  This 
new  array  is  distributed  in  nature.  The  improved  design 
not  only  reduces  the  internal  loss,  but  also  increases  the 
available  power  density  since  the  total  length  of  the 
array  is  reduced  dramatically.  The  performance  of  these 
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distributed  arrays  is  studied  here  both  by  computer 
simulation  and  measurement  on  actual  devices. 

II.  SIMULATION 

To  optimize  the  design  of  the  distributed  array,  w'e 
need  to  determine  Uie  junction  position  along  the 
microstrip  so  that  all  the  junctions  are  strongly  phase 
locked  with  a  nearly  constant  phase  with  respect  to  the 
RF  current  in  the  microstrip.  A  computer  simulation 
has  been  carried  out  to  study  the  behavior  of  such  an 
array.  The  simulation  is  based  on  the  perturbation 
theory[5]  of  the  Resistively  Shunted  Junction  (RSJ).  For 
a  A'  junction  one  dimensional  array,  if  the  junctions  arc 
dc  biased  in  series,  then  the  voltage  across  the  mtr,  junc¬ 
tion  can  be  expressed  as[o,6] 

'*’0  ' 

where  is  the  total  RF  current  going  through  the  mth 
junction  due  to  all  the  junctions  in  the  arrays 
.V 

11  +  C.M 

t  =  1 

and  wJJ,.  and  are  the  voltage,  the  down  con¬ 

version  coefficient  and  the  dynamical  resistance  of  the 
mth  junction  respiectively  in  the  absence  of  interactions. 

is  the  rf  voltage  amplitude  generated  by  the  Fth  junc 
tion.  The  interactions  among  junctions  are  included  in 
the  equation  through  the  complex  mutual  admittance 
with  phase  The  linear  component  of  the 

phase  of  each  junction  is  given  by  For  given 

ofierating  voltage  =  \'/I ^Ry  junction  parameters  -J;,. 
Q  and  (,  can  be  calciilated[5]. 

For  an  A’  junction  array,  Fxjs.  1  and  2  give  a  set  of 
2,V  coupled  equations.  To  study  the  phase  locking 
behavior,  these  2  A'  equations  must  be  solved 
simultaneously.  Phase  locking  is  achieved  once  the 
voltages  (9^  across  all  the  junctions  become  equal.  This 
locking  process  can  be  quite  complex  in  the  type  of 
resonant  structure  studied  here  since  Y can  varv 
rapidly  with  frequency  near  the  operating  fwint.  The 
structure,  as  shoiA’n  in  Fig.  la,  is  typical  of  those  studied. 
The  total  length  of  the  array  is  400 /jm  long,  about  one 
wave  length  at  an  operating  frequency  /  “  300  GHz 
The  20  junctions  in  the  array  are  separated  into  two 
groups.  The  inter-junction  spacing  within  each  group  is 
10  ^m,  and  the  inter-group  spacing  is  200 //m,  or  a  half 
wave  length  at  the  designed  frequency.  The  arrav  is 
tX)C  1993  IEEE 
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terminated  at  the  right  end  with  a  load  Z i  =■  30  0.  In 
the  computer  simulation,  identical  RSJs  are  assumed 
with  normal  resistance  =  0.5  0,  and  I^R^  ~  600 /i\'. 
Ihe  simulation  prtKess  starts  with  randomly  given  initial 
phases  for  all  the  junctions.  The  iteration  proceeds  until 
the  phase  locking  condition  is  achieved.  The  simulation 
results  reveal  several  important  features  of  distributed 
arrays. 

A.  Disiribuftons  of  rf  current  and  junction  phase 

The  distributions  of  RF  current  and  the  phase 
difference  A6  =  ^;«nc(ion  ~  l>etwMn  the  junction  and 
RF  current  at  320GHz  (/=1.18A)  are  shown  in  Fig.  lb 
when  all  the  junctions  are  phsise  locked.  The  horizontal 
axis  in  Fig.  lb  represents  the  locations  of  junctions  in  the 
array.  As  expected,  a  strong  resonance  is  established. 
The  magnitude  of  the  rf  current  follows  the  profile  of  a 
standing  wave.  It  can  be  noted  that  there  is  an 
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Fig.  1  (a)  The  structure  of  the  array,  (b)  .Vmplitude  of 
RF  current  (o)  and  the  phase  difference  Aff  (  V  ). 

appreciable  RF  current  (Jrj—  203 pA)  going  through  the 
load,  i.f.  a  traveling  wave  component  is  also  present 
near  the  load. 

B.  Stability  of  phase  locking 

The  stability  of  phase  hxking  is  usually  measured  by 
the  range  in  bias  current  over  which  phase  locking  can  be 
sustained[5].  It  has  been  shown  that  this  stability  of 
phase  locking  is  related  to  the  phase  difTerence  A5 
between  the  junction  and  the  RF  rurrent[5].  Once  phase 
locked,  the  phase  difference  is  always  within  0  and  n 
with  the  maximum  stability  at  ir/2  and  the  least 


stability  at  0  and  rr.  By  applying  this  stability  rule,  one 
can  observe  from  Fig.  lb  that  the  junctions  located  near 
the  RF  current  anti-nodes  are  the  most  stable.  The 
shaded  areas  indicate  the  unstable  regions  where  no  June, 
tions  should  be  placed  The  existence  of  these  less  stable 
regions  is  due  to  the  low  RF  current  amplitude  and  the 
less  ideal  phase  differences  both  of  which  .'educe  tie 
loc  king  strength.  • 

III.  INTEGRATED  SUB.MILLIMETER 
WAVE  RECEIVER 

The  above  simulation  results  outline  the  beisic  design 
rules  for  the  distributed  array.  Arrays  with  the  same 
parameters  as  those  shown  in  Fig.  la  have  been 
tabricated.  To  evaluate  their  performance,  we  performed 
both  RF  power  and  radiation  linewidth  measurements. 
In  order  to  measure  the  radiation  linewidth  in  the 
submillimeter  wave  range,  we  designed  a  prototype 
Integrated  Submillimeter  Wave  Receiver  (ISWR) 
schematically  shown  iu  Fig.  2.  Here  two  arrays,  a 
Josephson  effect  detector  and  an  SIS  mi.ver  are 
integrated  on  a  single  silicon  substrate.  One  array  func¬ 
tions  as  the  local  oscillator  and  pumps  the  SIS  junction. 
The  other  array,  weakly  coupled,  serves  as  the  signal 
source.  Both  arrays  have  the  same  design  parameters. 
The  load  for  these  arrays  is  the  SIS  mixer  with  a  normal 
state  resistance  /?„  between  20  to  40  Q.  As  sIkomi  m 


Fig.  2  Schematic  of  Integrated  Submillimeter  Wa^e 
Receiver. 

Fig.  2,  the  SIS  mixer  is  grounded  via  a  lower  impi  cUnte 
(  ~0.5n)  RSJ  junction.  The  primary  function  of  thn- 
single  RSJ  is  to  measure  the  rf  current,  therefore  ilie  rf 
power,  delivered  to  the  SIS  junction  from  the  local 
oscillator  array.  Also  included  on  the  chip  are  .m  IF 
filter,  a  coupler/dc  block  and  a  dc  bias  network  Tl.e 
characteristic  impedance  of  the  microstrip  transm,"'inn 
line  is  20  Q.  The  IF  signal,  coming  from  micri  -trip. 
passes  through  the  on-chip  IF  filter  and  goes  to  a  >0  il 
SMA  connector  via  a  50  G  coplanar  wave  guide  tr.Misi- 
tion.  These  impedances  arc  chosen  so  that  there  >  a 
reasonable  impedance  match  between  different  sec:  r,ns 
However,  no  effort  has  been  made  to  optimize  them. 
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In  ordff  that  the  two  arrays  ran  be  frequency  tuned 
ifldfpendcntly.  phase  locking  should  not  exist  between 
ihrm.  One  function  of  the  small  coupling  capacitor  in 
ihf  microstrip  is  to  isolate  the  two  arrays  so  that  the 
(nrasurrd  locking  range  corresponds  to  only  about  10% 
ofih'ir  frequency  difference  at  the  10  GHz  IF  frequency. 
The  required  LO  power,  based  on  the  optimum  pumping 
poner  equation  from  Tucker's  theoryjT],  is 

=  (  h«r,/e)^/2/?„  (3) 

where  -  is  the  LO  pumping  frequency  and  /?,  is  the 
normal  state  resistance  of  the  SIS  junction.  The 
dimensionless  variable  a  =  el'^f)/h«;,  is  a  measure  of  the 
LO  power  level.  Assuming  all  the  junctions  are  in  phase, 
one  can  estimate  the  power  delivered  to  the  SIS  mixer  as 


where  <  is  the  RF  voltage  generated  by  each  junction. 
This  implies  that  a  minimum  of  6  junctions  are  needed 
for  the  llOnW  required  pumping  power,  assuming 
=  0.5 n,  =  25n  and  a  ~  2.  The  value  of  a  is 
chosen  based  on  the  assumption  that,  the  conversion  is 
on  the  peak  of  the  first  Bessel  function  7,(o)(7]. 


The  fabricated  ISWR  is  designed  to  operate  around 
300  GHz.  The  RSJ  junction  area  is  2  x  2  pm^,  which 
requires  a  current  density  —  30  kA/cm^  for  the 
desired  /^  =  1  mA.  Low  induct.ance  (£.  CrO.lpH)  Au 
rcsi'iors  are  used  for  resistive  shunts.  For  the  SIS  mixer, 
the  product  is  designed  close  to  3.  With.  R„  =  20 
0  :0  0.  a  =  4  kA/cm^  is  necessary  for  the  I  x  1pm* 
junction  used.  Because  these  current  densities  are  quite 
different,  two  separate  Nb/AlO/Nb  tri-layers  with 
different  oxidation  times  were  used  The  detailed  .Nb 
junction  fabrication  process  can  be  found  in  reference  4. 
In  ihe  aesign,  contrary  to  the  simulation,  the  junctions 
in  the  array  are  parallel  biased  in  order  to  compensate 
for  the  non-uniformity  of  junction  parameters.  Each  10 
junction  section  of  the  array  can  be  bisised 
independently. 


IV.  EXPERIMENTS  AND  RESULTS 


All  the  experiments  wer:  done  in  a  p-metal  and  lead 
can  shielded  dewar  with  thi?  sample  immersed  directly  in 
liquid  helium.  All  the  junctions/arrays  were  powered  by 
batteries  with  careful  filtering.  The  output  power  from 
the  20  junction  LO  array  was  measured,  using  the  SIS 
junction  biased  at  8  mV,  where  it  has  a  resistance  of 
26  fi,  as  the  load.  The  Shapiro  step  width  of  the  single 
RSJ  detector  junction  in  series  with  the  SIS  junction  was 
used  to  determine  the  RF  current  through  this  load. 
Figures  shows  the  RF  power  measured  at  1.57  K  as  a 
function  of  LO  frequency.  The  mciximum  power  of  1.8 
pW  is  obtained  at  the  'iesigned  resonant  frequency  of 


300  GHz.  The  dash  line  is  the  expected  power  calculated 
by  assuming  the  in-phase  operation  of  all  the  junctions. 
Eq.4.  The  effect  of  junction’s  shunt  inductance  on  the 
RF  power  measurement  is  estimated  to  be  les.s  than  13% 
at  300 GHz.  The  data  show  that  the  LO  array  delivered 
full  power  to  the  load,  implying  that  all  the  junctions  in 
the  distributed  array  oscillate  coherently  and  in  phase. 

The  radiation  linewidth  is  measured  by  using  a  double 
heterodyne  mixing  technique.  The  first  mixing  is  done 
on  chip  using  the  ISVVR.  The  first  IF  output,  typically 
10  to  18  GHz,  is  then  brought  out  of  the  He  dewar  via 
semi-rigid  coaxial  cable.  It  is  mixed  with  a  second 
microwave  LO  giving  the  second  IF  at  100  MHz,  which 
is  coupled  to  the  deteetbr  through  an  IF  bandpass  filter 
whose  bandwidth  varies  from  1  MHz  to  4  MHz.  The  I-V 
curves  of  the  SIS  mixer  junction  are  shown  in  Fig.  4. 
Curve  a  is  the  unpumped  SIS  junction.  The  lack  of 


Fig.  3  Measured  RF  power  spiectrum  (o)  from  20  junc¬ 
tion  LO  array  delivered  to  SIS  junction  biased  at 
8  mV  (R^=26n)  at  1.57  K.  Dashed  line  is  the 
estimated  power  from  E^u.  4. 


Fig.  4  I-V  curves  of  SIS  with  10  junction  LO  array  off 
(a)  and  on  (b,  c). 

apparent  supercurrent  in  these  curves  is  probably  due  to 
the  trapped  magnetic  flux.  Traces  b  and  c  are  the 
pumped  1-V  curves  when  the  local  oscillator  array  was 
operating  at  210  GHz  and  312  GHz  respectively.  The 
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qu&si-parlicle  steps  ran  be  clearly  seen  both  below  and 
above  the  gap  voltage.  When  the  ISWR  is  operating  in 
the  receiver  mode  for  the  linewidth  measurement,  the  20 
junction  LO  array  delivers  too  much  power  to  the  SIS 
mixer.  To  ci't  down  the  pumping  ptower,  only  10  junc¬ 
tions  are  used.  K  seems  that  the  power  coupled  to  the 
SIS  mixer  at  210  GHz  is  still  higher  than  the  optimum 
level. 

During  linewidth  measurements,  the  signal  array  is 
biased  at  256  GHz  whereas  the  LO  array  is  biased 
at240GHz.  The  SIS  mixer  is  biased  at  the  center  of  the 
first  quasi-particle  step  where  the  mscximum  response  is 
expjected.  The  sjoectrum  of  the  first  IF  taken  at  1.57  K  is 


shown  in  Fig.  5.  At  this  temperature,  the  measured 
linewidth  —  10.2  MHz,  which  is  the  convolution  of 
the  linewidths  fro;r;  LO  and  SIG  arrays,  is  about  509? 
larger  than  the  linewidth  estimated  by  assuming  that  the 
noise  is  entirely  due  to  the  Nyquist  noise  of  shunt 
resistors.  At  4.0  K,  the  measured  linewidth,  23.1  MHz.  is 
15%  larger  than  that  estimated.  Both  the  low  frequency 
and  the  down  converted  high  frequency  noises  are 
included  in  the  estimate.  A  more  detailed  linewidth 
analysis  will  be  published  elsewhere. 

VII.  CONCLUSION 

A  single  chip  receiver  and  signal  source  operating  near 
250  GHz.  using  Josephson  efTecl  array  oscillators  and  an 
SIS  mixer  has  bezii  fabricated  and  tested.  The  IF. 
coupled  off  the  chip  to  a  room  temperature  spectrum 
analyzer,  permitted  a  direct  measurement  of  the 
oscillator  linewidth.  The  arrays  were  of  a  new  design, 
with  junction  distributed  over  a  substantial  part  of  a 
wavelength.  The  measured  jxjwer  and  linewidth  were  in 
good  agreement  with  those  obtained  from  computer 
simulations,  implying  complete  in  phase  locking  of  the 
sources. 
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The  radiation  linewidth  of  a  phase-locked  distnbuted  Josephson  junction  array  source  at  230 
GHz  has  been  measured  as  a  function  of  temperature  using  an  integrated  submillimeter  wave 
receiver/source  The  measured  total  linewidth  (including  signal  and  local  oscillator  arrays, 
having  10  junctions  each)  is  10  MHz  at  1.6  K.  The  measured  data  agree,  within  50% 
uncertainty,  with  the  linewidth  obtained  from  computer  simulation  assuming  the  Nyquist  noise 
current  spectral  density  generated  by  the  junction  shunt  resistors. 


The  development  of  sources  based  on  the  phase-locked 
arrays  of  Josephson  junctions'  ’  has  reached  a  sute  where 
a  useful  power  level’  (several  jiW)  can  be  delivered  to  a 
high  impedance  load  (20-600)  in  the  submillimeter  wave 
range.  However,  measurements  of  the  linewidth  for  these 
arrays  have  not  been  made  Previous  measurements  of  the 
radiation  linewidth  of  phase-locked  Josephson  junction  ar¬ 
rays  were  performed  between  2  and  20  GHz'*  on  the 
lumped  arrays,  w  here  the  length,  /,  of  the  region  containing 
junctions  was  much  smaller  than  the  wavelength.  A  The 
theoretically  predicted  dependencies  of  the  radiation  line- 
width  on  the  junction  dynamic  resistance  and  on  the  num¬ 
ber  of  junctions  were  o’oserved.  The  potentially  much  more 
useful  array  sources  operating  in  or  near  the  submillimeter 
wave  range,  however,  have  a  distnbuted  structure  (/>>!). 
In  this  letter,  we  report  and  compare  simulations  and  mea¬ 
surements  of  the  radiation  linewidth  of  phase-locked  dis¬ 
tnbuted  Josephson  junction  arrays  operating  at  250  GHz 
The  radiation  linewidth  of  a  single  resistively  shunted 
Josephson  junction  (RSJ)  has  been  well  studied,’  ’  and  the 
predictions  were  venfied  expenmentally  *  For  a  RSJ.  the 
noise  IS  dominated  by  the  Nyquist  noise  current  generated 
b)  the  shunt  resistor  R,  producing  a  linewidth 


ISw  =  IT 


(5^/(0)  +  la^S i((i>xi) ). 


(1) 


where  Rj  is  the  dynamic  resistance  of  junction  at  the  op¬ 
erating  point,  a  is  the  down  conversion  coelficient,  ’Pomh/ 
2e  is  the  flux  quantum,  and  5/(0)  and  Si(uq)  are  the 
current  noise  spectral  densities  at  low  frequencies  and  at 
the  Josephson  frequency  Wq.  respiectively,  and  are  given  by 
the  Nyquist  formula  including  zero  point  fluctuations^ 


5/(01) 


fka 

itR 


(2) 


For  an  array  with  ,V  junctions,  as  shown  in  Fig  1(a), 
there  are  ,V  noise  sources  each  of  which  contributes  to  the 
linewidth.  Therefore  the  total  linewidth  for  a  phase-locked 
array,  assuming  these  noise  sources  are  incoherent,  can  be 
written  as 


Aw=irj^j  I^/jLl5/«(0)  +  2ai^/„(a>o)] 

=  rr(^J  (5/(0)-(-2a’5,(£Uo))  (3) 

where  the  last  equal  sign  holds  for  identical  junctions.  Here 
R^i,s  {<P(/2ir)  (dioo/d/ i)  is  the  dynamic  resistance  of  the 
Jrth  junction  and,  in  general,  is  affected  by  the  interactions 
of  the  Jrth  junction  with  other  junctions  in  the  array  and 
with  the  circuit,  e  g.,  by  the  rf  current  7^/ going  through  the 
Jrth  junction.  One  of  the  great  advantages  of  phase-locked 
anays  is  that  the  oscillation  linewidth  is  considerably  nar¬ 
rowed  compared  to  the  linewidth  of  a  single  junction  os¬ 
cillator  identical  to  those  in  the  array  For  a  lumped  array, 
the  linewidth  can  be  reduced  by  a  factor  of  N  since  the 
junction's  dynamic  resistance  R^i,  can  be  1/<V  of  that  of  a 
single  junction.*  In  a  distributed  array,  however,  th.s  sim¬ 
ple  relation  may  not  hold  since  I^can  vary  with  position  in 
the  array,  so  must  in  general  be  determined  by  com¬ 
puter  simulation. 

We  have  numericaily  simulated  the  distribution  of  the 
junctions'  dynamic  resistances  for  the  distributed  array 
structures  shown  in  Figs.  1(a)  and  1(b).  In  these  arrays, 
20  junctions,  in  two  groups  of  10.  are  coupled  together 
through  a  20  n  microstnp,  which  is  designed  to  be  one 
wavelength  A  long  at  the  intended  operating  frequency  of 
300  GHz.  As  seen  in  Fig  1(b),  the  left  end  of  the  trans- 
missiim  line  is  open,  and  the  other  end  is  either  terminated 
with  a  load  (upper  array),  which,  in  this  case,  is  a 
superconductor-insulator-superconductor  (SIS)  mixer 
junction  with  a  normal  state  resistance  of  about  30  D.  or  is 
weakly  coupled  to  the  mixer  (lower  array)  through  a  ca¬ 
pacitor.  The  interjunction  spacing  is  10  jim,  and  the  center 
to  center  intergroup  spacing  is  200  jim.  i.e.,  A/2.  The  RSJs 
used  in  the  arrays  have  resistive  shunts  with  a  designed 
7^  =  600  jiV  and  7f  =  0  9  fl  With  these  parameters,  the 
normalized  voltage,  y/I^R,  is  close  to  unity  at  the  designed 
operating  frequency.  This  assures  that  the  radiation  from 
the  RSJ  IS  almost  monochromatic,  and  at  the  same  time 
the  phase  locking  among  junctions  is  the  strongest. 

The  computer  simulation  used  to  determine  R^t,  is 
based  on  a  perturbation  theory'  where  the  phase-locking  rf 
current  through  a  junction  in  the  microstrip  is  treated  as  a 
perturbation  to  the  analytical  solution  of  the  RSJ  moucl. 
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FIG  1  (a)  Structure  of  disiribuied  lotephson  junction  array  dmgncd  lo 
have  a  resonance  at  }00  GHz  The  right  end  is  terminated  with  a  load  of 
)0  n  (bl  Layout  of  integrated  submillimeter  wave  receiver/source  The 
signal  array  (bottom)  is  weakly  coupled  to  the  SIS  mixer  The  /F  is 
coupled  off  the  chip  and  measured  by  a  room-temperature  spectrum  an- 
al>zer  For  clanty.  only  the  dc  bias  of  one  section  of  the  parallel  biased 
SIG  array  is  shown  (dashed  line). 

tahich  is  denoted  through  superscript  u.  eg,  ui“(I).  The  rf 
current  is  related  to  the  junctions'  oscillations  through 
the  frequency  dependent  admittances^  The  .V 

phase-locking  equations  for  the  frequencies  of  the  ;unvtion 
oscillation  can  be  written  as 

)  =‘uJ(/4  )  -t- 

Xcos()-*,„-r-d„-0i),  (4) 

where  e*  is  the  magnitude  of  the  fundamental  Josephson 
oscillation,  and  ;S  the  linear  component  of  the  kth  junc¬ 
tion's  phase  In  order  to  calculate  the  djnamic  resistance 
of  indoidual  junctions  in  the  phase-locked  array,  the 
dc  voltage  change  was  determined  as  the  bias  current 
through  the  Ath  junction  was  varied  while  holding  all  other 
biases  constant  In  Fig  2,  the  dynamic  resistances 
(solid  squares),  together  with  the  rf  cunent  (open  cir¬ 
cles).  are  plotted  as  a  function  of  junction  position  in  the 
array  operating  at  a  frequency  such  that  the  wavelength 
A=0  %l  As  can  be  seen,  the  dynamic  resistance  varies  from 
junction  to  junction,  reflecting  the  distributed  nature  of 
this  array  The  average  of  RJ|,^i^  mfl  is  slightly  less  than 
that  of  R'ii/S  =  34  mO  expected  from  a  lumped  array  due 
to  the  resonant  structure  of  the  arrsy  Also  shown  in  Fig. 
2  are  the  calculated  Rji  for  both  the  upper  (down  solid 
triangles)  and  the  lower  arrays  (up  solid  triangles)  when 
only  the  group  of  10  junctions  nearer  the  mixer  are  biased. 

To  measure  the  radiation  linewidth  in  the  submiliime- 
ter  wave  range,  we  have  developed  an  integrated  submilli¬ 
meter  w ave  rcceiver/source  ( ISWR/S ) .  ’  This  contains  two 
arrays,  having  the  same  dimensions  as  shown  in  Fig  1(a). 
and  one  SIS  mixer  junction,  fabricated  using  Nb/A10,/Nb 
technology  The  upper  array  serves  as  the  local  oscillator 
(LO)  and  pumps  the  SIS  mixer,  the  other  (weakly  cou¬ 
pled)  serves  as  the  signal  source  (SIG)  All  of  the  junc¬ 
tions  were  powered  by  low  noise  battery  operated  power 
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FIG  2  Calculxtfd  junction  dynamic  rcsistancr  u  a  function  of  junction 
position  in  the  LO  array  when  20  junctions  were  used  (solid  squares ).  and 
in  the  LO  array  (down  solid  tnangles).  and  in  the  SIG  array  (up  solid 
inan|les)  with  only  the  10  junctions  nearer  the  mixer  operating  The  rf 
current  (open  circles),  when  20  junctions  are  used,  is  also  shown 

supplies,  which  gave  a  negligible  contribution  to  the  line- 
width  While  the  diffierent  groups  of  junctions  could  be 
biased  independently,  junctions  within  a  group  were  biased 
in  parallel  from  a  common  source  through  a  senes  of  in¬ 
terlocking  dc  superconducting  quantum  interference  de¬ 
vice  loops'  ’  having  a  high  impedance  at  oiq.  All  voltage 
leads  were  filtered  by  an  8-ms  RC  filter,  whereas  the  high 
current  bias  leads  were  filtered  by  14  ms  RC  filter. 

The  power  spectral  density  of  the  radiation  was  mea¬ 
sured  by  using  a  double  heterodyne  mixing  technique  The 
signal  from  the  SIG  array,  operating  at  256  GHz.  was 
mixed  down  to  an  //■=  16  GHz  by  the  on  chip  £*S  mixer, 
which  was  pumped  by  the  LO  array  oscillating  at  240 
GHz  One  reason  for  selecting  this  relatively  high  IF  fre¬ 
quency  was  to  ensure  that  the  interaction  of  the  two 
weakly  coupled  arrays  did  not  affect  their  linewidth  To 
obtain  the  proper  power  level  from  the  LO  array,  only  the 
groups  of  10  junctions  closer  to  the  SIS  m’xer  were  actually 
used  The  IF  was  coupled  via  coaxial  ct.ble  out  of  the  cry¬ 
ostat  where  it  was  further  mixed  dowii  to  110  MHz  and 
detected.  An  isolatoi,  inserted  Jt  the  top  ol  the  cryostat, 
prevented  signals  and  noise  from  the  room-temperature 
electronics  from  coupling  back  to  the  sample.  The  band¬ 
width  of  the  second  //^bandpass  filter,  1  MHz  for  T <  2  5 
K,  and  4  MHz  for  T>2.5  K,  determined  the  frequency 
resolution  in  the  spectrum  measurements  The  detected 
power  has,  in  the  worst  case,  a  peak  signal  to  background 
noise  ratio  larger  than  5.  This  background  noise  was  sub¬ 
tracted  from  the  measured  spectrum.  The  solid  circles  in 
Fig  3  show  the  measured  total  linewidth  (LO  array  plus 
SIG  array)  as  a  function  of  liquid  helium  bath  temper-iiure 
from  4  0  to  1.5  K. 

In  order  to  accurately  calculate  the  expected  linewidth. 
one  needs  to  know  the  dynamic  resistance  of  each  junciion 
as  well  as  the  current  noise  spectrum.  Our  parallel  bias 
scheme  made  it  impossible  to  directly  measure  the  clyn.imu: 
resistance  of  individual  junctions  However,  the  dynamic 
resistances  of  the  arrays  could  be  measured  and  used  to 
csnmaie  the  junctions'  dynamic  resistances  based  on  the 
simulation  The  dynamic  resistance  of  the  l.O  (SIG)  array 
varies  from  65  mfl  ( 58  mfl )  at  4  2  K.  to  55  mil  ( 40  mSl  I 
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rlG  3  Measured  radiation  linewidth  (aolid  ctrriet)  aa  ■  function  of 
temperature  The  solid  line  is  the  calculated  total  linewidth  The  long 
dashed  line  is  the  coninbution  from  low  fre<iuenc>  current  noise  The 
*hor:  dashed  line  ts  the  contribution  from  high  frequenc>  current  noise 
including  tero  poini  huctuations  at  Josephson  oscillation  frequency 


at  !,57  K.  The  linewidths  of  the  individual  arrays  are  then 
calculated  from  £q  (3)  using  the  calculated  values  of 
shown  III  Fig.  2(b),  scaled  to  give  an  average  equal  to  the 
measured  parallel  dynam.c  resistance.  This  introduces  an 
estimated  uncertainty  of  about  30^c  for  the  calculated  line- 
width.  An  additional  uncertainty  in  the  companson  of  the 
calculated  and  measured  linewidth  is  introduced  by  the 
cor  volution  of  the  LO  and  SIG  linewidths  dunng  the  on 
r.nip  mixing.  The  finite  linewidth  of  the  LO  array  affects 
the  total  linewidth  through  the  change  of  mixer  gam. 
which  is  mainly  a  function  of  the  LO  power,  within  the 
linewidth  of  the  LO  array.  Around  the  usual  applied  LO 
pumping  power  level  ( <  the  optimum  LO  pumping 
power),’  the  mixer  conversion  gain  is  roughly  proportional 
to  the  applied  LO  power. '“Such  a  relation  gives  a  convo¬ 
luted  linewidth  equal  to  the  sum  of  the  two  linewidths. 
assuming  Lorcntzian  line  shapes  The  nonlinear  gam- 
power  relation  of  the  physical  mixer  could  introduce  an 
error  of  about  \0^c  m  the  convoluted  linewidth  The  solid 


line  in  Fig.  3  shows  total  linewidth  calculated  from  Eqs. 
(2)  and  (3)  using  these  approximations. 

Our  measured  linewidths  from  10  junction  distributed 
arrays  sources  agree  with  those  expected  from  the  Nyquist 
noise  due  to  the  junctions'  shunts  within  the  50%  uncer¬ 
tainty  of  out  companson  This  indicates  that  additional 
noise  sources  such  as  the  shot  noise  of  tunnel  currents  or 
coherent  rf  noise  currents  coupled  via  the  coaxial  cable  are 
not  dominate.  The  analysis  presented  can  therefore  serve  as 
a  guide  to  the  expected  linewidths  in  the  much  larger  ar¬ 
rays  being  developed  as  practical  submillimeter  sources 
where,  for  constant  array  impedance,  one  expects  Aoict  1/ 
.Vf 
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COMPLETE  PHASE-LOCKING  IN  A  ONE-DIMENSIONAL  SERIES 
BIASED  JOSEPHSON-JUNCTION  ARRAY 

Siyu&n  H&n,  A.  H.  Worah&m,  and  J.  E.  Lukens 
Department  of  Physica 
SUNY  at  Stony  Brook 
Stony  Brook,  NY  11794 


Abstract  ~Tht  rf  power  at  —  300  GHz  coupled  to  a 
52  Cl  load  by  a  one-dimensionzil  dc  seriea  biased 
Josephson  junction  array  has  been  measured  using  an  on- 
chip  Josephson  detector.  From  the  measured  rf  power 
and  dc  voltage  levels  we  conclude  that  the  100  junctions 
in  this  ID  array  were  locked  in-phase.  The  effects  of  the 
finite  inductance,  associated  with  the  junction’s  shunt 
resistor,  on  array’s  output  power  and  detector’s  current- 
voltage  characteristics  are  also  discussed. 

I.  INTRODUCTION 

Pha»^  locked  Josephson  junction  arrays  have  many 
advantages  over  single  junction  oscillators  as  millimeter 
and  submillimeter  sources.  Compared  to  a  single 
junction  of  shunt  resistance  and  characteristic  voltage 

s  I^R,  where  is  critical  current,  a  one-dimensional 
(ID)  rf  series  array  of  total  resistance  R^  consisting  of  N 
junctions,  each  having  the  same  characteristic  voltage 
as  that  of  the  single  junction,  is  expected  to  increase  the 
rf  power  delivered  to  a  matched  load  Ri  =  R,  by 
times,  and  to  reduce  the  radiation  linewidth  by  a  factor 
of  while  preserving  the  rapid  tuning  speed  of  the 
single  junction  source  [1-3].  All  of  these  properties  of  the 
ID  array  have  been  experimentally  demonstrated 
recently  [4-7]  in  arrays  having  several  tens  of  junctions. 
For  exsunple,  the  rf  power  delivered  to  a  60  fl  load  at 
350  GHz  was  observed  to  be  equal  to  the  maximum 
expected  power  of  a  completely  phase  locked  array  [4]. 
(Recently  published  results  of  computer  simulations  of 
ID  arrays  [8]  showing  a  maximum  power  far  less  than 
the  coherent  sum  of  the  powers  from  the  individual 
junctions  are  not  applicable  to  our  work,  since  these 
simulations  did  not  include  any  phase-locking  inter¬ 
actions  among  the  junctions.) 

II.  EXPERIMENT 
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Fig.  1.  Schematic  of  the  100-junction  series  biased  ID 
array  (x  represents  the  Josephson  junction). 


since  it  requires  much  lower  bias  currents  [9]  and  permits 
the  use  of  various  shunting  schemes  for  line-width 
reduction.  Here,  we  report  new  results  on  a  ID  array  of 
senes  ic  biased  junctions  which  was  fabricated  at  IBM 
by  the  Consortium  for  Superconducting  Electronics 
(CSE)  using  the  Planarized  All-Refractory  Technology 
for  low-Tg  Superconductivity  (PARTS)  developed  by 
Ketchen  et.  al  [lOj.  This  new  technology  produced  a 
100-junction  ID  array  with  a  spread  in  of  about 
±3%. 


For  the  experimental  results  referred  to  above,  the 
scatter  in  the  junction  parameters  within  the  arrays 
made  it  necessary  to  use  a  parallel  dc  bias  scheme  to 
achieve  phase  locking  in  these  rf  series  arrays  [4,5].  Our 
simulations  show  that,  for  example,  a  scatter  in  of 
less  than  ±  7%  is  required  to  achieve  phase  locking  with 
series  dc  biasing.  Scries  biasing  is  never-the-lcss  desirable 


This  100-junction  ID  array  was  designed  to  test  both 
the  long  standing  predictions  for  phase-locking  with 
seriea  bias  and  whether  the  current  fabrication 
technology  for  Nb  based  Josephson  junctions  could 
provide  adequate  junction  uniformity  for  complete  phase 
locking.  As  seen  in  the  schematic  of  the  sttray  in  Fig.  1, 
the  junctions  are  placed  at  intervals  of  Aq  =  370  pm 
along  a  20  0  serpentine  microstrip  transmission  line. 
Both  ends  of  the  transmission  line  are  terminated  by 
quarter  wave  (Ajj/4)  stubs.  The  detector  junctii.n  has 
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the  same  nominal  critical  current  and  shunt  resistance  as 
the  junctions  in  the  array.  For  such  an  array  of  ideal 
resistively  shunted  junctions  operated  at  friquency 
/ss  Wp/Ag,  where  is  the  phase  velocity  of  the 

transmission  linu,  the  maximurri  ptower  delivered  to  a 
resistive  load  is 


^  2(R,  +  rJ 


(1) 


where  R^^  NR]  if  ail  the  junctions  are  identical.  Eq.  1 
is  exauit  only  for  average  dc  bias  voltages  of  each 
junction  V  ~ 

The  power  delivered  to  the  load  is  measured  by 
the  amplitude  of  the  constant  voltage  step  in  the 
detector’s  dc  current-voltage  characteristics  (IVC).  The 
critical  current  is  measured  from  the  IVCs  of  the  array 
to  be  80  at  ~  1.4  K.  The  total  resistance  of  the 
array  R,  =  162  is  measured  from  the  IVC  with  the 
Josephson  current  completely  suppressed  by  applying  a 
n  ignetic  field,  giving  the  mean  shunt  resistance  per 
junction  Rj  —  1.62  n.  The  damping  parameter 

0^  S  2‘irV ^R jCI^Q  of  the  junctions  is  estimated  to  be 
~  0  07,  where  C  is  the  junction  capacitance  and 
#0  s  A/fe  is  the  flux  quantum.  The  IVC/junction  of 
the  100  junctions  in  series  and  the  IVC  of  the  detector 
are  shown  in  Fig.  2.  These  IVCs  are  nearly  identical 
showing  the  very  good  uniformity  in  and  R j  among 
all  the  junctions  (  ~  3%).  In  Fig.  3  the  detector’s  IVC, 
for  the  array  biased  at  62.4  mV,  is  shown  to  clearly 


I  (mA) 


Fig.  2.  IVCs  of  the  array  and  the  detector.  The  artay’s 
voltage  is  shown  as  the  total  voltage  divided  by 
the  number  of  junctions  (A^=100).  For  clarity, 
the  detector’s  IVC  (dttshed  line)  is  shifted 
vertically. 


display  a  current  step  at  624  ±0.2  >iV,  t.e.  at 
amplitude  of  the  rf  current  derived 


from  the  size  of  the  current  step  (n=l  Shapiro  step)  is 
I^j  —  93  /iA.  The  value  of  =  52  fl.  as  determined 
from  the  IVC  of  the  load  resistor,  thus  gives  the  rf  power 
delivered  to  the  load  resistor  as  0  22  /iW. 
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Fig.  3.  Detector  junction's  IVC  with  array  biased  at 
V'j,rroy=02.4  mV  (solid  line)  and  array  off 
(dashM  line). 


The  value  of  s  93  txk  is  derived  using  a  modified 
Resistively  Shunted  Junction  (RSJ)  model-  Resi.stively- 
Inductively  and  Capacitively  Shunted  Junction  (RICSJ) 
model  -  for  the  detector  junction  which  takes  into 
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Fig.  4  (a)  Schematics  of  a  Resistively-Inductively 

Capacitively  Shunted  Junction  (RICSJ)  and  tb) 
its  equivalent  circuit. 


1  1 


2491 


account  the  finite  inductance  associated  with  the  external 
shunt  resistor  placed  across  the  tunnel  junction. 
According  to  RSJ  model  the  size  of  the  first  Shapiro  step 
(  ~  34  fjiA)  requires  the  detector  Junction  be  irradiated 
an  rf  current  having  200  /iA  amplitude,  giving  the  rf 
power  coupled  to  the  load  resistor  to  be  about  1  ijW. 
Bowever.  the  existence  of  the  [xarasitic  inductance 
associated  with  the  shunt  resistor  (cf  Fig.  4) 
substantially  increases  the  fraction  of  the  external  rf 
current  flowing  through  the  tunnel  junction  when  the 
frequency  f>(RJL,)/2r.  This,  in  turn,  increases  the 
step  amplitude  for  a  given  rf  amplitude  (in  the  linear 
response  range).  Fig.  5  shows  the  relationship  between 


Islep/Ie 


Fig.  5  Size  of  nssl  Shapiro  step  /,,ep  ''*■  appl'et^  rf 
current  amplitude  7,^  of  the  RSu  model  (open 
triangles)  with  /#q  =  4.8V^i.,  J^=0.07  and  the 
RICSJ  model  (open  circles)  for  =  0.22  with 
same  /  and  where  =  2x1, /^/^g 

1 and  for  the  two  different  junction  models  based 

J  JICp  ** 

on  our  numerical  simulations.  A  more  detailed  analysis 
of  the  RICSJ  model  will  be  published  elsewhere. 

III.  DISCUSSION 


where  a  =  2'rfL^/Rj.  This  increase  in  the  rf  power 
coupled  to  the  load  is  not  hard  to  understand  since  the 
inductor  provides  considerable  impedance  at  high 
frequencies  (2x/>  7£,/I.,)  hence  reducing  the  rf  piower 
dissipated  by  the  junction’s  shunt  resistor.  In  deriving 
Eq.  2  we  have  made  the  assumption  that  tunnel 
junaion’s  quasi-particle  conductance  <  ^/Rj-  for 
the  tunnel  junctions  in  the  ID  array  this  is  justified. 
According  to  Eq.  2  the  rf  power  coupled  to  R^  at 
~  300  GHz  should  increase  to  0.12  ^W,  still  much  lower 
than  the  observed  value. 

One  of  the  possible  explanations  is  that  when  rf 
current  is  close  to  or  greater  than  the  critical  current  of 
the  junctions  and  the  structure  is  resonant  type,  ^he 
perturbation  approach  and  the  lumped  circuit 
approximation,  which  leads  to  Eqs.  1  and  2,  is  not 
adequate.  It  has  been  experimental)  •  observed  and 
confirmed  by  the  numerical  simulation  that  in  a  ID 
array  the  rf  current  through  'he  array  can  be  several 
times  /j.  due  to  the  strong  coupling  between  series 
junctions  embedded  in  a  microstrip  line  and  the 
resonator  structure  of  the  ID  array  [1].  This  explanation 
is  supported  by  the  resonant  structure  observed  in  the 
ID  array's  IVC  (cf.  Fig.  6). 
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This  measured  rf  power  coupled  to  the  load  is  more 
than  twice  the  value  of  P /ij  =  9b  nW  calculated  from 
Eq.  1,  using  independently  measured  RSJ  parameters. 
However,  Elq.  1  is  derived  under  the  assumption  that  the 
junctions  shunt  inductance  is  zero.  In  our  sample, 
each  junction  has  an  inductance  ~  0.9  pH  in  series 
with  the  shunt  resistor  (cf  Fig.  4).  Taking  into  account 
this  inductance,  the  rf  power  coupled  to  the  load  would 
be  increased  to 


1  +[<>/(  WR^/fl,)P 


(2) 


Fig.  6  Part  of  the  IVC  of  the  array  showing  the 
resonant  structures. 


While  the  maximum  measured  rf  power  is  about  twice 
that  predicted  (in  the  approximation  used)  we  have  a 
high  degree  of  confidence  in  this  result.  In  principle, 
independent  of  the  junction  and  array  model,  the  total 
power  dissipated  by  all  the  junctions  and  the  load  at  dc 
and  the  rf  frequency  should  not  exceed  the  input  power 
supplied  by  the  dc  current  bias  source.  As  a  check,  we 
carefully  measured  the  input  power  and  the  dc  power 
dissipation  and  always  find  sufficient  excess  power  to 
supply  the  rf  power  to  tlie  load. 
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IV.  SUMMARY 

In  summary,  we  have  demonstrated  that  complete 
phase  locking  has  been  achieved  in  a  senes  biased  ID 
array  fabricated  using  PARTS  technology.  While  the 
design  was  restricted  to  rather  low  rf  power  due  to  the 
relatively  low  critical  current  density  (  ~  3  kA/cm^) 
which  was  required  for  the  rest  of  the  devices  (e.g.  dc 
SQUIDs)  on  the  same  wafer,  the  measured  rf  power 
coupled  to  the  52  Q  load  is  0.22  pW,  exceeding  the 
maximum  power  expected  from  the  array  by  pertur¬ 
bation  analysis.  This  observed  rf  fxjwer  delivered  to  the 
load,  along  with  the  excellent  agreement  between  the 
bias  voltage  of  the  array  and  the  frequency  of  the 
detected  radiation,  provide  strong  evidence  for  complete 
in-phase  locking  of  the  junctions  in  this  100  junctions  ID 
array  with  series  dc  bias. 
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As  millimeter  and  submiLUmeter  sources,  phtuse-locked  array  of  N 
Josephson  junctions  can  overcome  many  shortcomings  of  the  single  junction 
source  such  as  low  source  impedance,  low  output  power,  and  large  linewidth.  In 
recent  years  it  has  been  demonstrated  that  small  linear  arrays  {N  <  100)  are 
capable  of  delivering  severed  fiW  power  to  a  load  of  10  to  60  fl  in  the  frequency 
range  of  ~  200  GHz  to  ~  300  GHz  with  linewidth  Af/f  <  5  x  10'®  [1-4].  For 
many  applications  one  would  like  to  have  higher  power  and  narrower  linewidth 
which  can  be  achieved  by  increasing  the  number  of  jimctions  deployed  in  the 
array  since  the  output  power  increase  as  and  linewidth  decreases  as  l/jV^  for 
arrays  having  NRj  =  Ri,  where  Rj  and  Ri  are  the  junction’s  shunt  resistance  and 
the  load  resistcince,  respectively.  There  are  two  basic  designs  for  the  medium  size 
{N ~  10^)  and  large  size  (iV~  10^)  linear  Josephson  array:  the  qucisi-lumped  array 
and  distributed  array  [4,  5].  In  a  quasi-lumped  array  the  disteince  between  two 
adjacent  junctions  is  Aq  (wavelength  at  designed  operating  frequency).  In  a 
distributed  array  a  group  of  m  jimctions  are  closely  packed  together  and  the 
disteince  between  two  adjacent  groups  is  Ag.  In  both  types  of  array  the  junctions 
are  embedded  in  the  treinsmission  line  structure  which  provides  the  long-range 
high-frequency  electromagnetic  couplings  between  the  junctions.  This  coupling  is 
essential  to  achieve  mutual  phcise  locking  in  the  linear  array.  Compared  to  the 
quasi-lumped  array,  the  distributed  array  requires  less  space  and  has  smaller 
internid  rf  loss  (especially  as  /g  approaches  the  superconducting  gap  frequency  of 
~  700  GHz  for  Nb). 

The  output  rf  power  coupled  to  a  load  R^  from  an  linear  array  of  N 
junctions  can  be  written  as 


,  __  hNVcfR, 

^  2{NRj  +  Rl)^ 


(1) 
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when  the  source  impedance  matches  to  the  load  Eq.  (1)  becomes 


where  =  Ic^j  is  the  characteristic  frequency  of  the  junction,  Iq  is  the  critical 
current  of  the  junction.  7  depends  on  the  dc  bias  voltage  and  is  usually  less  than 
unity. 

We  have  designed  medium  size  (1^=500)  distributed  array  (m  =  10) 
oscillators  for  operation  around  400  GHz.  Each  resistively  shunted  Nb/AlO/Nb 
tunnel  junction  is  5.2  /im*  in  size  and  has  designed  parameters  of  /c=1.4  mA 
{Jc  =  27000  A/cm^)  eind  Rj=  0..5  fl  which  give  =  0.7  mV, 

2;  0.25  and  =  2TrIcLg/^Q  ~  0.8,  where  Cj  and  Lg  are  the  shunt  capacitiince  of 
the  junction  and  the  parasitic  inductance  associated  with  the  shunt  resistor, 
respectively.  $q  is  the  flux  quantum.  The  total  length  of  the  array  is  about  50  Ap 
(at  ~  400  GHz).  The  Nb  microstrip  line  is  250  nm  thick  and  8  fim  wide  with  670 
nm  thermally  evaporated  SiO  between  the  microstrip  line  cind  the  340  nm  thick 
Nb  ground  plane.  This  gives  Zq  c;  12  Cl.  The  designed  output  power  to  the  50  Cl 
load  is  about  100  /iW  around  400  GHz. 

The  Nb/AlO/Nb  trilayer  films  were  made  at  AT&T  Bell  Laboratory  and 
the  arrays  were  fabricated  at  IBM  using  the  PARTS  process  developed  by 
Ketchen  et.  al  [6].  The  measured  parameters  of  a  typical  single  junction  from 
wafer  are  /<;  =  2.5  mA,  Ry  =  0.56  11  and  mV  at  4.2  K.  This 

corresponding  to  a  Jc  of  about  50000  A/cm  which  is  about  twice  of  the  designed 
value.  This  results  in  a  value  of  /3c  —  0.56  and  1.6  which  are  quite  far  away 
from  the  designed  value  and  unfavorable  for  phase-locking.  The  measured  load 
resistance  is  62  Cl. 

The  output  rf  power  to  the  load  resistor  has  been  measured  using  the  on- 
chip  Josephson  junction  detectors.  The  cimplitude  of  the  rf  current  through  the 
detector  junction  is  estimated  by  '■ompanng  the  size  of  the  first  Shapiro  step  and 
the  suppression  of  the  critical  current  to  those  obtained  from  the  numerical 
simulations  using  measured  Vq,  Ic,  and  R j.  The  estimated  output  rf  pv/wer  P,y  is 
45  ±  15  /iW  using  the  suppression  of  and  is  50  ±  25  using  the  size  of  the 
first  Shapiro  step.  The  measured  rf  power  is  consistent  with  =  65  /iW 
calculated  from  Eq.  (1)  [7].  The  measurement  of  the  output  power  spectrum  Pr/^f) 
is  currently  in  progress. 
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In  summary,  a  linear  Josephson  effect  array  oscillator  with  series  dc  bias 
has  been  designed  to  deliver  rf  power  of  about  100  /rW  to  a  50  fl  load.  The 
measured  output  power  to  the  62  Cl  load  at  ~  400  GHz  is  50  /iW,  lower  than 
designed  value.,  due  to  the  difference  between  the  target  critical  current  and  that 
of  the  sample. 
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For  many  submillimcter  wave  (SMW)  frequency  applications  it  is  essential  to  have  a  compact 
source  capable  of  delivering  rf  power  in  the  range  of  10—  100  /iW  to  a  load  around  50  f2.  The  nighest  rf 
ptower  measured  from  a  Josephson  effect  oscillator  was  7  /iW  at  ~  ,300  GHz  reported  by  Wan  el.  ai  [1- 
4].  Recently,  we  have  obtained  rf  power  of  greater  than  50  /iW  coupled  to  ~  CO  G  load  around  400  GHz 
from  distributed  linear  Josephson  arrays  [4]  of  500  junctions  with  series  dc  biasing. 

The  rf  power  coupled  to  a  load  Rj^  from  an  completely  phase-locked  linear  array  of  N  Josephson 
junctions  can  be  written  as 

P^  =  K(^NVfRJ[2{NRj  +  R^)'^]  (1) 

Here.  Rj  is  the  normal  state  resistance  of  each  junction,  =  I^R j  is  the  characteristic  voltage  of  the 
junction  and  I ^  is  the  critical  current.  7  and  k  are  factors  of  order  unity  which  depend  on  the  junction 
parameters  and  array  design.  In  general  7  and  «c  must  be  obtained  by  computer  simulation. 

The  distributed  array  (shown  in  Fig.  1)  has  500  resistively  shunted  junctions  placed  along  a 
microstrip  transmission  line.  Ten  junction.«  are  grouped  into  a  lump  in  which  two  adjacent  junctions  are 
separated  by  10  /rm.  .Adjacent  lumps  are  separated  by  one  wavelength  at  the  designed  operating 
frequency.  Load  resistors  and  Josephson  junction  detectors  are  placed  at  both  ends  of  the  array.  The 
arrays  were  fabricated  at  IBM  using  P.4RTS  technology  [5]  With  high  critical  current  density 
(,/j  SI  30  -  50  kA/cm^)  .Nb/AIO,/Nb  trilayers  made  at  AT&T  Bell  Laboratory.  700  —  800  nm  of 
thermally  evaporated  SiO  was  used  as  the  dielectric  between  the  Nb  microstrip  and  the  300  nm  ,Nb 
ground  plane.  Each  load  resistor  and  detector  junction  has  independent  dc  bias  leads  which  allows  us  to 
make  four  terminal  measurement  to  determine  their  parameters.  Rj  is  measured  by  suppressing  the 
junction's  critical  current  with  a  magnetic  fielo.  The  capacitance  of  the  6  junctions  is  estimated  to 
be  ~  290  fF  including  parasitic  capacitance.  The  parasitic  inductance  L  associated  with  the  external 
shunt  resistor  is  0.18  pH  estimated  from  the  physical  dimensions  of  the  shunt. 

The  coherent  rf  current  generated  by  the  array  is  measured  by  the  detector  junctions  placed 
after  the  loads.  The  magnitude  of  the  n=I  Shapiro  step  [6j  is  used  to  obtain  the  amplitude  of  the 
coherent  rf  current  through  the  junction  and  the  load  bv  fitting  the  measured  IV  curve  to  that  of 
numerical  simulations  using  measured  junction  parameters.  For  numerical  simulations,  the  detector 
junc'.icn  has  to  be  described  by  the  Resistively-lnductively  and  Capacitively  Shunted  Junction  (RICSJ) 
model  due  to  the  rather  large  value  of  LI^  [7J. 

Five  samples  have  been  completed  and  tested  for  SMW  generation.  AH  of  them  have  >  20  /iW 
coupled  to  a  60  D  load  around  400  GHz.  The  dc  IV  curves  of  a  500  junction  array  and  its  detector 
junction  from  sample  W11D8  are  shown  in  Fig.  2.  The  array  has  an  average  critical  current  as  3.2  mA 
and  an  average  Rj  as  0.4  fl.  The  detector  junction  has  as  3.45  mA,  Rj  as  0.4  Q,  Fig.  3  shows  the 
measured  IV  curve  around  the  n=:l  Shapiro  step  at  1.9  K  and  that  of  the  simulation  using  '.he  measured 
parameters  above  and  =  1.22  mA.  The  measured  load  resistance  is  68  Q  for  this  sample  thus  giving 
rf  powe.'  ~  50  pW  at  400  GHz.  The  size  of  the  n=l  Shapiro  step  of  sample  W8F5  as  a  function 

of  array's  oscillation  frequency  is  displayed  in  Fig.  4.  Note  that  the  array  oscillator  generated  usable  rf 
power  (  >  1  /iW)  into  ~  60  fl  over  a  fairly  wide  frequency  range  (300  -450  GHz).  .More  than  10  pVV  rf 
power  at  500  GHz  has  been  obtained  from  sample  WIID8. 

In  summary,  series  dc  biased  500  junction  distributed  Josephson  junction  arrays  have  been 
designed,  fabricated  and  tested.  The  measured  maximum  coherent  rf  power  coupled  tc»  a  60  fl  load  at 
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400  GHz  is  greater  than  50  /iV\’.  Coherent  rf  power  above  10  pW  has  also  been  measured  at  500  GHz 
from  the  array  oscillator. 
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Fig.  1  Schematics  of  500  junction  senes  biased 
distributed  array. 


Fig.  2  dc  rV’  curves  of  the  array  and  the  de¬ 
tector  junction.  Array's  voltage  has 
been  divided  by  500  (sample  WT1D8). 


Fig-  3  The  data  and  simulation  of  detector’s  IV 
curve  with  array  on  (sample  W11D8). 


Fig.  4  Measured  size  of  n=:l  Shapiro  step 
vs.  frequency  (sample  W8F'5). 
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ABSTRACT 

A  submillimeter  wave  source  based  on  a  new  design  using  Josephson  junction 
arrays  has  been  developed  and  tested.  The  maximum  rf  power,  delivered  to  a  68  f?  load 
and  detected  on-chip,  was  47  pW  at  394  GHz.  .Significant  power  was  detected  at  a 
number  of  frequencies  from  300  GHz  to  500  GHz  where  the  power  was  10  pW.  The 
observed  power  at  the  designed  operating  frequency  near  400  GHz  is  consistent  with  all 
500  junctions  in  the  series  biased  array  delivering  current  in-phase  to  the  loads.  This  is 
in  agreement  with  simulations  of  smaller  arrays  of  the  same  design.  The  linewidth. 
inferred  from  the  measured  resistance  at  the  point  of  maximum  power,  with  T  —  4.2  K.  is 
less  than  1  MHz.  The  minimum  inferred  linewidth  near  400  GHz,  at  somewhat  lower 
power,  is  about  100  kHz. 
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Many  submilliineter  wa’’?  (S.MW)  frequency  applications,  such  as  radio  astronomy 
and  high  speed  communication,  require  a  compact  SMW  source  capable  of  delivering  rf 
poi».cr  in  the  range  of  lO-lOOpW.  For  frequencies  above  300  GHz,  the  lack  of  compact 
sources  is  particularly  serious  requiring  increasingly  inefficient  multiplication  of 
millimeter  wa”e  sources.  Lately  there  has  been  a  substantial  effort  to  develop  Josephson 
effect  sources,  which  can  work  to  at  least  1  THz  [1],  for  applications  in  the  SMW  [2-4], 
The  most  promising  results,  recently  obtained,  have  come  from  three  types  of  sources: 
long  single  junctions  flux  flow  oscillators  [2],  various  types  of  two  dimensional  arrays  of 
small  junctions  [3]  and  one  dimensional  (ID)  small  junction  arrays  [4],  discussed  here.  In 
genereJ  a  power  of  about  1  pW  is  obtained  from  these  sources,  although  in  many  cases  the 
source  impedance  is  less  than  1  fl.  The  highest  rf  power  previously  reported  from  a 
Josephson  effect  oscillator  has  been  7  ftW  at  ~  300  GHz  by  Wan  et.  al  [4]  using  a  ID 
array.  Here  we  report  results  from  an  improved  design  for  Josephson  SMW  sources — 
distributed  linear  arrays  of  500  Josephson  junctions  with  series  dc  biasing — generating  an 
rf  power  of  47  pW  euound  394  GHz  and  10  pW  at  500  GHz  into  an  on-chip  load  of  6S  Cl. 

Small  resistively  shunted  junctions  (RSJ),  described  by  the  RSJ  model  [5],  are 
characterized  by  a  critical  current  1^  and  a  shunt  resistance  Rj,  and  generate  a  peak  rf 
voltage  (V^^l^Rj)  with  a  source  impedance  Z,  at  the  Josephson  frequency. 

/;  =  4S3.6  GHz/m\’.  For  bias  voltage  *1  and  Z,--*Rj.  The  junctions  are  placed 

in  a  superconducting  microstrip  and  interact  via  the  rf  current  I,j  generated  in  the 
microstrip  by  the  Josephson  oscillations.  For  properly  designed  arrays  [6]  the  rf  voltages 
across  the  junctions  can  be  made  to  add  in  series.  This  will  be  true  independent  of  the 
relative  polarities,  along  the  microstrip,  of  the  junctions’  biases,  i.  e.  the  dc  voltage  can 
sum  to  (series  bias)  or  to  zero  (parallel  bias).  The  power  from  such  an  array  (see 
Fig.  1 )  is  then  just, 


2{NRj  +  Rl?  ■ 


(1) 


Here  N  is  the  number  of  rf  series  elements  in  the  array,  and  R^  is  the  load  impedance 
seen  by  the  array,  k  =  1  is  the  ideal  case,  but  is  in  general  somewhat  less  than  unity  due 
to  losses  in  the  coupling  structure  or  imperfect  phase  alignment  of  the  junctions.  For  ID 
arrays  (Fig.  la-c),  N  is  just  the  number  of  junctions  in  the  array.  For  2D  turays  (Fig. 
Id)  of  the  types  so  far  reported,  the  rf  current  flow  is  ideally  along  rows  of  junctions. 
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with  the  phase  of  the  rf  current  transverse  to  this  flow  (i.e.  across  a  column)  being 
constant.  Equation  1  can  then  be  applied  to  2D  arrays  as  well  if  one  replaces  N  by  N,. 
the  number  of  series  elements  (i.e.  columns)  in  the  array  and  Rj  by  Rj,  the  resistance  of 
cl  series  element.  The  power  from  an  array  matched  to  the  load,  i.e.  with  =  Ri,  is 

Pn  ~  {^,Vc)^/BRi,  and  can  in  principle  be  increased  just  by  incretising  N,  while  keeping 
NRj  and  V).  constant.  This  implies  that  a  N,  and  ~  Rild/S. 

Increasing  the  source  power,  therefore,  requires  increasing  7^,,  adding  together 
enough  series  elements  to  match  the  array  impedance  to  the  load  and  maintaining  phase 
coherence  throughout  the  array.  In  practice,  there  are  two  essential  limitations  to  this 
increa.se.  First,  when  1^,  becomes  greater  than  several  milliamps,  there  exists  the 
potential  for  flux  flow  instabilities  which  would  dissipate  power  internally  in  both  single 
junctions  and  in  the  columns  of  2D  arrays.  It  is  not  yet  clear  whether  there  is  an 
advantage  to  ID  or  2D  arrays  for  maximizing  /<.,  while  maintaining  phase  stability.  The 
second  problem,  common  to  both  ID  and  2D  arrays,  is  that  when  the  length  of  the  array 
becomes  an  appreciable  fraction  of  the  wavelength,  care  must  be  taken  to  maintain  the 
proper  phase  relationship  between  the  junction  oscillations  and  the  rf  current  throughout 
the  array.  This  is  required  so  that  all  junctions  phase  lock  in  a  stable  manner  and 
provide  current  to  the  load  with  the  same  phase.  Especially  for  the  SMW  source,  this 
problem  must  be  solved  to  permit  the  use  of  enough  junctions  to  obtain  power  levels 
much  above  a  microwatt.  It  is  this  second  problem  which  we  address  in  this  paper. 

Fully  coherent  ID  arrays  have  been  demonstrated  using  both  a  lumped  circuit 
design  (Fig.  la),  where  the  array  length  L  <  A/10  [6]  (A  is  the  wavelength  in  the 
microstrip  coupling  the  junctions),  or  the  so-called  quasi-lumped  design  where  the 
junctions'  spacing  is  A  [4].  The  number  of  junctions  in  both  designs  is  severely  limited  in 
the  SMW.  A  maximum  of  about  N  =  10  can  be  used  in  the  SMW  for  lumped  arrays 
made  with  standard  fabrication  technologies.  For  quasi-lumped  arrays,  losses  in  the 
microstrip  coupling  the  junctions  as  well  as  the  larger  array  size  are  limitations.  Much 
denser  packing  can  be  achieved  if  the  junctions  can  be  distributed  along  the  microstrip 
for  a  significant  fraction  of  A.  Figure  Ic  shows  such  a  structure  where  groups  of  M 
junctions  are  placed  at  intervals  of  A  or  A/2  along  the  microstrip.  Our  computer 
simulation  shows  that  for  same  N  and  M  arrays  with  A/2  inter-lump  distance  have  a 
wider  tuning  range. 

Distributed  arrays  (shown  in  Fig.  2)  of  500  junctions  have  been  designed. 
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fabricated,  and  tested  for  SM\N’  generation.  Resistively  shunted  Nb/AlO^c/Nb  tunnel 
junctions  are  placed  in  groups  of  ten  junctions  {M  —  10)  with  adjacent  junctions  in  the 
group  separated  by  10  fim.  Adjacent  groups  are  separated  by  one  wavelength  at  the 
designed  primary  operating  frequency  of  ~  390  GHz.  Load  resistors  and  Josephson 
junction  detectors  to  measure  the  rf  currents  are  placed  at  both  ends  of  the  array.  The 
arrays  were  fabricated  at  IBM  using  Plemarized  All-Refractory  Technology  for 
Superconductivity  [7]  with  high  critictJ  current  density  (J^ss40-60  kA/cm^) 
Nb/A10,/Nb  trilayers  made  at  ATicT  Bell  Laboratory.  700  nm  of  thermcdly  evaporated 
SiO  was  used  as  the  dielectric  between  the  Nb  microstrip  tind  the  300  nm  Nb  ground 
pl^me,  which  was  placed  on  top  of  the  array.  Each  load  resistor  and  detector  junction  has 
independent  dc  bias  leads  allowing  their  parameters  to  be  determined  using  four  terminal 
measurements.  Rj  is  measured  by  suppressing  the  junctions’  critical  current  with  a 
magnetic  field.  The  capacitance  of  the  6  junctions,  including  parasitic  capacitance, 
is  estimated  from  the  geometry  and  the  specific  capacitance  of  46  fF/^m^  to  be  ~  290  fF 
.  The  parasitic  inductance  L  associated  with  the  external  shunt  resistor  is  0.18  pH 
estimated  from  the  physical  dimensions  of  the  shunt.  The  junctions  in  the  array  are 
biased  in  series  using  a  common  dc  current.  This  requires  a  higher  degree  of  junction 
uniformity,  but  is  important  for  larger  arrays.  For  example,  the  arrays  studied  here 
would  have  required  a  bias  current  of  over  2  A  for  parallel  bias  instead  of  the  5m.A 
actually  used. 

The  coherent  rf  current  generated  by  the  array  is  measured  by  the  detector 
junctions  placed  after  the  lumped  resistor  loads.  The  magnitude  of  the  n  =  l  Shapiio  step 
[8]  is  used  to  obtain  the  amplitude  of  the  coherent  rf  current  I,j  through  the  detector 
junction  and  the  load  resistor  by  fitting  the  measured  1\’  curve  to  that  of  numerical 
simulations  using  measured  junction  parameters.  For  numerical  simulations,  the  detector 
junction  has  to  be  described  by  the  Resistively  inductively  and  Capacitively  Shunted 
Junction  model  [9]  due  to  the  rather  large  value  of  LI^. 

Five  samples  have  been  completed  and  tested  for  SMW  generation.  All  delivered 
more  than  20  fiW  of  rf  power  near  394  GHz  to  loads  of  about  65  fh  For  sample  W'llDS, 
the  a^•erage  critical  rurrent  of  the  junctions  in  the  array  is  7^  =  3.2  mA  and  the  averag*' 
junction  resistance  is  Rj  -  O.SSfh  The  detector  junctions  have  an  7^  ~  3.5  niA,  and  a 
7?j  =  0.38  0.  The  maximum  power  from  the  array  occurs  at  394  GHz,  however, 
significant  power  is  av  ailable  at  a  number  of  frequencies  in  the  300  GHz  to  500  GHz  range 


with  10 ^\V  of  power  delivered  to  the  loads  at  500GHz.  Figure  3  shows  the  measuied  /V' 
curve  of  a  detector  around  the  n=l  Shapiro  step  with  the  array  operating  at  394  GHz  and 
T  =  1.9K.  This  is  comp2ired  with  that  of  the  detector  simulation  using  the  measured 
junction  parameters  and  an  rf  current  amplitude  from  the  ^ray  of  =  1.18mA.  The 
measured  total  load  resistance  is  68  fl  for  this  sample  giving  an  rf  power, 
~  47 /iW  at  394  GHz.  Since  the  total  array  resistance  is  190(1  (the  actual 
junctions  shunt  resistors  were  somewhat  greater  than  the  design  vadues)  the  maximum 
avciilable  power  from  this  array  (into  a  190  fl  load)  would  be  64  pW  at  394  GHz. 

The  measured  power  of  ~50//\V  at  394  GHz  compares  reasonably  with  the  value 
of  96  pW  from  Eq.  1,  assuming  perfect  alignment  of  the  junctions'  phases  ztnd  no 
transmission  line  losses  (i.e.  k  =  1).  However,  k  is  in  general  always  less  than  unity  due 
to  imperfect  phase  alignment  and  the  losses  in  the  coupling  structure  (in  our  case,  a 
microstrip  transmission  line).  The  measured  value  of  k  =  47/96  ~  0.49  at  394  GHz  can 
be  entirely  accounted  for  b}'  a  surface  resistance  of  7  mfl,  assuming  nearly  perfect  phase 
alignment.  O-  computer  simulations  of  an  array  with  A'  =  100  but  otherwise  identical 
in  structure  tr,  that  measured  indeed  show  a  negligible  phase  spread  among  the  junctions. 
Although  the  surface  resistance  of  the  actual  sample  has  not  been  measured,  the  value 
inferred  from  the  rf  power  measurement  (7  mfl)  is  consistent  with  the  measured  surface 
resistance  at  about  400  GHz  in  N'b/SiO/Xb  microstrip  resonators  [10]  and  compares 
fa%orably  to  that  obtained  by  Cucolo  ft  al.  [llj. 

The  second  major  advantage  of  the  phase-locked  array  source  is  that  the  radiation 
linewidfh  can  be  sub.^tantially  reduced  since  the  phase-locking  suppresses  the  frequency 
modulation  of  the  array  due  to  the  random  noise  sources  associated  with  the  junction 
shunt  resistors.  The  linewidth  of  a  series  array  of  A’  identical  Josephson  junctions  is 
approximately  given  by  [0.  12,  13]: 

“  (^ )  2Q'xcoth(r')]  ,  X  -  ^  (2) 

where  kg  is  the  Boltzmann  constant,  h  is  Planck’s  constant,  4'q  is  the  magnetic  flux- 
quantum,  Rj  is  the  dynamic  resistance  of  an  independent  junction,  and  a  ~  0.42  is  the 
down  conversion  coefficient.  Equation  2  requires  some  modification  in  distributed  arrays, 
since  the  strength  of  the  locking  can  depend  on  the  position  of  the  junction.  However, 
direct  measurements  of  the  linewidth  [l3]  in  a  small  array  consisting  of  two  groujis  of 


r 


junctions  separated  by  A/2  show  that,  within  a  factor  of  2  to  3,  Eq.  2  provides  a  reliable 
estimate  of  Af.  Using  vedues  of  Rj  and  Rj  as  the  average  over  the  array,  the  implied 
linewidth  at  394  GHz  is  about  730  kHz  at  7*  =  4.2  K  decreasing  to  475  kHz  at  1.6  K,  where 
I  ~  5.9,  so  zero  point  fluctuations  dominate.  The  value  of  Rj  varies  with  bias  current 
giving  a  minimum  implied  linewidth  of  about  100  kHz  at  T  =  4.2K  near  375  GHz  but 
with  a  somewhat  lower  power  level. 

In  summary,  500  junction  Josephson  effect  arrays  using  dc  series  bias  and  a 
distributed  junction  layout  have  been  designed,  fabricated  2md  tested.  The  measured 
maximum  coherent  rf  power  coupled  to  a  68  fl  load  at  394  GHz  is  47  pW,  implying  an 
available  power  of  64p\V.  An  rf  power  of  10  pW  has  also  been  measured  at  500  GHz 
from  the  arraj’  oscillator.  At  4.2  K  the  linewidth  near  400  GHz,  calculated  from  Eq.  (2) 
using  the  measured  array  parameters,  is  less  than  730kHz. 

This  work  is  supported  in  part  by  BMDO-IST  through  the  RADC  and  by  OXR, 
with  support  for  the  sample  fabrication  provided  by  CSE,  which  is  supported  in  part  by 
ARPA.  We  gratefully  thank  R.E.  Miller  and  W.H.  Mallison  at  AT<k"T  for  providing  the 
Nb/AlO/Nb  trilayer  films  and  M.  Ketchen’s  group  at  IBM  for  fabricating  the  sample. 
W’e  also  thank  A.W’.  W'orsham  and  A.  Jain  for  their  valuable  assistance  and  discussion. 


35 


REFERENCES 


[1]  R.P.  Robertazzi  and  R.A.  Buhrman,  IEEE  Trans.  Mag.  MAG-25,  13S4  (19S9). 

[2]  \'.M.  Zhang,  D.  Winkler,  and  T.  Claeson.  Appl.  Phys.  Lett.  62,  3195  (1993).  \’.P. 
Koshelets,  A.\’,  Shchukin,  S.\'.  Shitov,  IEEE  Trans.  Appl.  Supercon.  3,  2524  (1993). 

[3]  S.P,  Benz  and  C.J.  Burroughs,  Appl.  Phys.  Lett.  58,  2162  (1991).  J.S.  Martens,  A. 
Pance,  K.  Char,  L.  Lee,  S.  Whiteley,  and  \'.M.  Hietala,  ibid.  63,  1681  (1993). 

[4]  K.  Wan,  A.K.  Jain,  and  J.E.  Lukens,  Appl.  Phys.  Lett.  54,  1805  (1989). 

[5]  D  E.  McCumber,  J,  Appl.  Phys.  39,  2113  (1968).  W.C.  Stewart,  Appl.  Phys.  Lett. 
12,  277  (1968). 

[6]  A.K.  Jain,  K  K.  Likharev,  J.E.  Lukens,  and  J.E.  Sauvageau,  Phys.  Rep.  109,  309 
(1984).  J.E.  Lukens,  in  Superconducting  Devices,  edited  by  S  T-  Ruggiero  and  D..A. 
Rudnian  (Academic.  San  Diego,  1991),  p.  135. 

[7]  M.B.  Ketchen,  D.  Pearson.  A.  Kleinsasser.  C.-K.  Hu.  M.  Smyth,  J.  Logan.  K. 
Stawiasz.  M  Jaso.  K.  Petrillo.  M.  Manny.  S.  Basavaiah.  S.  Brodsky.  S.B.  Kaplan. 
W.J.  Gallagher,  and  M.  Bhushan,  Appl  Phys.  Lett.  59.  2609  (1991).  R.E.  Miller, 
W  H.  Mallison.  A  W.  Kleinsasser.  K..A.  Delin.  and  E.M  Macedo.  Appl.  Phys.  Lott. 
63.  1423  (19931, 

[S]  S  Shapiro,  Phys.  Rev.  Lott.  11.  80  (1963). 

[9]  S,  Han,  .A  H  \\'oishair..  and  J.E  Lukens.  IEEE  Tran>  .Appl  Su]>ercon.  3.  2489 
(1993). 

[10]  B.  Bi.  K.  Wan.  W.  Zhang.  S  Han.  and  J.E  Lukons.  IEEE  Trans.  .Aj^pl.  Supercon 
1.  145  (1991). 

[11]  .A.N.  Cucolo,  S.  Pace,  R.  A’aglio.  \’.  Lacquaniti.  and  G  Marullo.  IEEE 
Trans.  .Magn,  MAG-17,  812  (1981). 

[12]  K.K.  Likharev  and  \'.K.  Semenov.  JETP  Lett.  15,  442  (1972). 

[13]  B.  Bi,  S.  Han,  and  J.E.  Lukons.  Appl,  Phys  Lett.  62.  2745  (1993). 


It) 


FIGURE  CAPTIONS 


Fig.  1.  Some  typical  structures  for  Josephson  junction  arrays:  (a)  ID  Lumped  array,  (b) 
Quasi-lumped  array,  (c)  Distributed  tu'ray.  (d)  2D  array.  A  is  the  wavelength  at  the 
primary  operating  frequency;  is  the  dc  bias  current. 

Fig.  2.  Schematics  of  the  A  =  500,  M  —  10  serial  biased  ID  distributed  Josephson 
junction  array  for  SMW  generation. 

Fig.  3.  The  measured  IV  curve  (solid  line)  ^l^ound  the  n  =  1  Shapiro  step  at  798  pV  of 
the  detector  junction  on  sample  \\11D8  under  394  GHz  irradiation  from  the  array  and 
the  71’  curve  obtained  from  computer  simulation  (solid  circles)  using  RICSJ  model.  The 
parameters  used  for  the  simulation  are  7^  =  3.47  mA,  C  =  294  fF,  I  =  0.178  pH. 
7?j  =  0.377  n,  and  7^^  =  1.18  mA.  This  rf  current  amplitude  corresponds  to  a  power  of 
47  p\V  in  the  load  of  68  H. 
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MISSION 

OF 

ROME  LABORA  TORY 


Mission.  The  mission  of  Rome  Laboratory  is  to  advance  the  science  and 
technologies  of  command,  control,  communications  and  intelligence  and  to 
transition  them  into  systems  to  meet  customer  needs.  To  achieve  this, 
Rome  Lab: 


a.  Conducts  vigorous  research,  development  and  test  programs  in  all 
applicable  technologies; 

b.  Transitions  technology  to  current  and  future  systems  to  improve 
operational  capability,  readiness,  and  supportability; 

c.  Provides  a  full  range  of  technical  support  to  Air  Force  Materiel 
Command  product  centers  and  other  Air  Force  organizations; 

d.  Promotes  transfer  of  technology  to  the  private  sector; 

e.  Maintains  leading  edge  technological  expertise  in  the  areas  of 
surveillance,  communications,  command  and  control,  intelligence,  reliability 
science,  electro-magnetic  technology,  photonics,  signal  processing,  and 
computational  science. 

The  thrust  areas  of  technical  competence  include:  Surveillance, 
Communications,  Command  and  Control,  Intelligence,  Signal  Processing, 
Computer  Science  and  Technology,  Electromagnetic  Technology, 
Photonics  and  Reliability  Sciences. 


